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STUDY OF THERMODYNAMIC VENT AND SCREEN 


BAFFLE INTEGRATION FOR ORBITAL STORAGE 
AND TRANSFER OF LIQUID HYDROGEN 
By E. C. Cady 

McDonnell Douglas Astronautics Company 


SUMMARY 

A comprehensive analytical and experimental program was performed to 
determine the feasibility and desirability of integrating an internal thermo- 
dynamic vent system (TVS) and a full wall-screen liner (WSL) for the orbital 
storage and transfer of liquid hydrogen (LH 2 ). First, the conceptual design 
of the tankage system was determined, then a comprehensive survey of 
screens for use in the WSL was performed. Ten screens, spanning a wide 
range of weaves and retention performance (bubble point), were selected for 
analytical and experimental study. The experimental study determined the 
screen bubble point, flow-through pressure loss, and screen roughness- 
induced pressure loss along rectangular channels lined with screen on one 
side, for the 10 selected screens, using LH 2 saturated at 34. 5 N/cm.2 
(50 psia). The screen bubble-point and flow-through pressure loss data 
agreed well with previously determined data. The channel pressure loss data 
were unique and were correlated with friction factor and Reynolds number 
using a roughness parameter based on screen wire diameter. 

The experimental data were used in an analysis to determine the flow and 
performance characteristics of the WSL annulus during low-gravity LH 2 out- 
flow. The TVS pump system was optimized for pumping characteristics and 
minimum weight. These analyses were combined to determine the optimum 
integrated system in terms of minimum weight. The optimum system char- 
acteristics were determined and system fluid-dynamic feasibility was 
established for six tanks ranging from 141. 6 m3 (5, 000 ft^) to 1. 416 m3 
(50 ft^) for orbital storage times of 30 and 300 days, and for several TVS and 
LH 2 transfer flowrates. 


1 




INTRODUCTION 


Future space systems will require feed systems capable of in-orbit 
storage, expulsion, and resupply of cryogens in a controllable and predictable 
manner. Various concepts have been proposed to accomplish in-orbit trans- 
fer of fluids using surface tension principles. Some of these concepts have 
been evaluated to a limited extent. For example, during small-scale experi- 
ments conducted on Apollo 14 (ref. 1), a flow model using a relatively simple 
perforated annular baffle successfully controlled fluid behavior during low- 
gravity transfer. However, to achieve similar success with liquid hydrogen 
(LH 2 ) during orbital storage and transfer, the heat transfer to the stored 
LH^ must be controlled. 

A number of techniques have been proposed to achieve the required 
thermal control. One concept uses a dual-screen liner, and is designed to 
hold the LH 2 off the tank wall to limit the heat transfer (ref. 2). This 
approach entails complexity in construction and relies on passive, g -dependent 
thermal control which has not been demonstrated in low gravity. Active sys- 
tems for thermal control, based on thermodynamic venting phase conversion, 
have been extensively developed under NASA contracts (refs. 2 and 3). These 
thermodynamic vent systems use a pump -mixer to obtain fluid -dynamic and 
heat-transfer processes that are not significantly g-dependent, and have been 
satisfactorily demonstrated in extensive ground tests. 

Proper integration of this proven low-gravity venting and thermal control 
system with a single -wall screen liner for liquid acquisition could provide a 
simply constructed, reliable, and proven solution to low-gravity LH 2 storage 
and outflow. Further, optimization of the thermodynamic vent system and 
single-wall screen liner configuration and flow characteristics could provide 
thermal and fluid dynamic control in the LH^ tank during inflow. Inflow control 
problems were studied by MDAC during Project THERMO (ref. 4) and were 
found to be critical for most in-orbit resupply systems. 

The overall system concept studied is shown schematically in Figure 1. 

The system consists of two major components: a single- screen complete wall 

liner and a pump- driven thermodynamic vent system. The annulus between 
the screen and the tank wall remains full of LH 2 at all times and serves two 
functions. First, it provides liquid communication from the outflow line to 
the bulk LH 2 in the tank, which, although its orientation in the tank is 
unknown, will certainly be in contact with the tank wall because of the wetting 
characteristics of LH 2 . This communication allows outflow and LH2 trans- 
fer in low gravity. Second, the annulus provides a flow path for pumped LH 2 
which will absorb tank incident heating, flow through the standpipe, and 
reject the absorbed heat to the thermodynamic vent system. 

To determine the important characteristics of this integrated thermo- 
dynamic vent and screen system, a three-task study was conducted. In 
Task 1, details of the system configuration were determined, all available 
screens were surveyed, and pertinent data compiled, and then 10 different 
screens were selected for use in the next two tasks. In Task 2, important 
characteristics of the 10 selected screens, such as bubble point, flow-through 
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loss, and channel-flow Joss were experimentally determined using LH-, 
saturated at 34. 5 N/cm (50 psia). In Task 3, the Task 2 experimental data 
were used in a theoretical analysis to determine the optimum system charac- 
teristics in terms of minimum weight for storage times of 30 and 300 days 
in orbit. 
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CONCEPTUAL DESIGN AND SCREEN SELECTION 


Tankage Design Characteristics 

Six basic tank sizes and configurations were evaluated throughout the 
study. The characteristics of these tanks are shown in Table I. The annulus 
gaps which give residual volumes of 1% to 5% of the tank volume were basic 
parameters used both in the Task 2 experimental work and in the Task 3 
analysis. Referring to Figure 1, details of the system are shown which are 
common to all of the tanks studied. The thermodynamic vent system (TVS) is 
mounted on a manhole to allow removal. A slipjoint in the central standpipe 
facilitates TVS removal and allows use of different materials for the standpipe 
(the tank and screen liner are assumed to be type 304 stainless steel). 

Because of the very low head requirements for the TVS pump, an axial flow 
type of pump, driven by an electric motor was chosen, and shown together 
with vanes to direct the TVS flow into the annulus gap. Solid baffles are used 
at each end of the standpipe to direct the TVS flow and the outflow or incoming 
flow from or into the annulus gap. The TVS vent flow is boiled in a heat 
exchanger tube bonded to the outside of the standpipe to minimize standpipe 
pressure drop. The TVS vent flow is used to cool the inflow line, and some 
of the TVS flow is directed through the inflow baffle to eliminate a hot spot at 
the inflow line. The detailed operation of the system is described in the sec- 
tion on Analytical Evaluation of the Tankage System. 


Screen Survey 

The screen used in the wall-screen liner (WSL) will have an important 
effect on the design of the integrated system. In surveying screen materials 
for possible application in the WSL, the important characteristics are screen 
weave, bubble point, and pressure drop for flow through the screen. The 
bubble point is defined as that pressure, or head, which can be supported by 
the vapor -liquid interface in the pores of the screen before vapor bubbles 
break through into the liquid, or "breakdown" occurs. The supported head, 

H, at a g -level, g, for cryogenic propellants (zero contact angle) has been 
shown to be: 


H = 


4>'cr 

gpD 


( 1 ) 


where <j>' = 4 for circular pores of diameter, D, (ref. 5). To properly 
evaluate equation (1), the fluid properties of surface tension, tr, and density, 
p, must be accurately specified. The latest data on LH£ surface tension from 
the National Bureau of Standards (NBS) (ref. 6) indicate a variation with tem- 
perature as shown in Figure 2. Accepted data on LH 2 density from the NBS 
(ref. 7) is shown in Figure 3. Equation (1) is plotted in Figure 4 for 
10. 1 N/cm^ (14.7 psia) and 34. 5 N/cm^ (50 psia) saturated LH 2 properties 
for three specified g-levels, 10 “1, 10"^, and 10“3 g's. Superimposed on the 
plot (extrapolated to 0. 1 g's for convenience) are LH 2 bubble-point test data 
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Volume, L/D Dia, Length, Area 

m^(ft' 5 ) m(ft) m(ft) rn (ft^) 


141.6 3.66 14.65 168.5 

( 5 , 000 ) 4 ( 12 . 02 ) ( 48 . 06 ) ( 1 , 814 ) 

14.16 1.70 6.80 36.3 

( 500 ) 4 ( 5 . 58 ) ( 22 . 31 ) ( 391 ) 

14.16 2.21 4.42 30.8 

( 500 ) 2 ( 7 . 25 ) ( 14 . 50 ) ( 331 ) 

14 . 16 3 . 00 3 . 00 28 . 3 

( 500 ) 1 ( 9 . 85 ) ( 9 . 85 ) ( 305 ) 

1.416 1.03 2.06 6.6 

( 50 ) 2 ( 3 . 37 ) ( 6 . 74 ) ( 71 . 2 ) 

1.416 1.39 1.39 6.1 

( 50 ) 1 ( 4 . 57 ) ( 4 . 57 ) ( 65 . 6 ) 


Gap width -cm (in. ) 
for residual of 


1 % 

2 % 

3 % 

4 % 

5 % 

0 . 84 
( 0 . 33 ) 

1 . 68 
( 0 . 66 ) 

2 . 52 
( 0 . 99 ) 

3 . 35 
( 1 . 32 ) 

4 . 19 
( 1 . 65 ) 

0 . 38 
( 0 . 15 ) 

0 . 79 
( 0 . 31 ) 

1 . 17 
( 0 . 46 ) 

1 . 55 
( 0 . 61 ) 

1 . 96 
( 0 . 77 ) 

0 . 46 
( 0 . 18 ) 

0 . 92 
( 0 . 36 ) 

1 . 37 
( 0 . 54 ) 

1 . 85 
( 0 . 73 ) 

2 . 31 
( 0 . 91 ) 

0 . 51 
( 0 . 20 ) 

0 . 99 
( 0 . 39 ) 

1 . 50 
( 0 . 59 ) 

2 . 01 
( 0 . 79 ) 

2 . 52 
( 0 . 99 ) 

0 . 213 
( 0 . 084 ) 

0 . 427 
( 0 . 168 ) 

0 . 643 
( 0 . 253 ) 

0 . 856 
( 0 . 337 ) 

1 . 069 
( 0 . 421 ) 

0 . 231 
( 0 . 091 ) 

0 . 465 
( 0 . 183 ) 

0 . 696 
( 0 . 274 ) 

0 . 930 
( 0 . 366 ) 

1 . 161 
( 0 . 457 ) 
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EQUIVALENT CIRCULAR PORE DIAMETER - MICRONS 


Figure 4. LH 2 Static Head Retention Versus Pore Diameter 




































(taken in 1 g) of various screens, as summarized in Table II. The only LH 2 
bubble-point test data known prior to the experimental work performed in this 
study were determined by MDAC (ref. 8) and Lockheed Missiles & Space Co, 
(LMSC) (ref. 9). The equivalent circular pore size, D, in microns which 
satisfies equation (1) with <t>' = 4 and the experimental H is shown in Table II 
for each screen. Also shown, are the screen manufacturer's absolute micron 
rating for each screen. For the Dutch weave material, the absolute rating is 
determined experimentally, usually by filtering a slurry of glass beads; the 
largest bead passed by the screen determines the absolute rating. For 
square-weave screens, the absolute rating is determined analytically as the 
largest inscribed circle in the weave pore. 

To determine the anticipated LH 2 retention height from the manufacturers 
absolute rating (in the absence of bubble point data) for a particular screen, 
the absolute rating pore diameter was used in equation (1) and cjd was deter- 
mined, based on the experimental LH 2 retention data. The equivalent cj>’ is 
close to 3 (rather than 4), with a deviation of generally less than 10% (see 
Table II and Figure 5). Therefore, for fine mesh screens in LH 2 , the bubble 
point can be expressed as 


H = 


3 cr 

g P P a 


( 2 ) 


where p a is the absolute rated pore diameter. 

The pressure drop for flow through a screen can be described in terms of 
friction factor, f, and Reynolds number, R, in the manner of Armour and 
Cannon (ref. 10). The correlation is 

f=-|+P (3) 


v AP € 2 D% 

where R = — B , f = — , a and (3 are experimentally determined 

pa 2 D QbpV 2 

constants (see symbols) and Q is a tortuosity factor (1. 0 for square 
weave screens, 1. 3 for Dutch weave screens), while viscosity, p, and 
density, P are fluid characteristics, and Y is the fluid approach velocity 
to the screen. 

For the laminar flow regime, where R is small, (3 is generally much 
smaller than a/R and can be ignored, resulting in f = substitution of f and 

R ix 

gives: 

AP e 2 Dg c „ 2 

pa D 

QbpV^ ' P V 
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TABLE II. - SCREEN BUBBLE POINT DATA 


Mesh 

1-G LH 2 

Source 

Experimental equivalent 

Manufacturer 

' s rated 


Deviation 


Bubble -Point 


circular pore size. 

absolute pore size. 


from <j>' = 3. 0, 


m (ft) 


(microns) 

(mic 

rons) 



% 

450 x 2, 750 

1. 151 (3. 775) 

MDAC 

9. 6 

7 

a. 

b, c* 

2. 92 

- 

2. 7 

325 x 2, 300 

0. 933 (3. 06) 

MDAC 

12. 0 

10 


b, c 

3. 33 

+ 

11. 0 


0. 884 (2. 90) 

MDAC 

12. 5 




3. 20 

+ 

6.. 7 

325 x 1,900 

>. 805 ( >2. 64) 

MDAC 

13. 7 


- 


- 



250 x 1,370 

>. 666 ( >2. 185) 

MDAC 

16.6 

12. 5 

a. 

b, c 

3. 01 

+ 

0. 3 

200 x 1, 400 

0. 582 (1. 91) 

MDAC 

19 

14 

a. 

c 

2. 95 


1. 7 





15 

b 


3. 16 

+ 

5. 3 

165 x 800 

0. 25 (0. 82) 

MDAC 

44. 5 

35 

b 


3. 15 

+ 

5. 0 





37 

d 


3. 32 

+ 

10. 7 

200 x 600 

0. 244 (0. 80) 

MDAC 

46 

30 

b. 

d 

2. 61 

_ 

13. 0 





40 

e 


3. 48 

+ 

16. 0 

80 x 700 

>. 247 (>0. 81) 

MDAC 

45 

35 

b 


3. 11 

+ 

3. 7 





40 

a. 

c 

3. 56 

+ 

18. 7 

325 x 325 

0. 186 (0. 61) 

MDAC 

60 

45 

b 


3. 0 


0 


0. 173 (0. 566) 

LMSC 

61. 5 

50. 8 

f 


3. 3 

+ 

10. 0 

50 x 250 

0. 123 (0. 405) 

LMSC 

86 

64 

f 


2. 98 

_ 

0. 7 


0. 112 (0. 367) 

MDAC 

100 

65 

d 


2. 60 

- 

13. 3 

200 x 200 

0. 117 (0. 384) 

LMSC 

90 

73. 6 

f 


3. 27 

+ 

9. 0 


0. 104 (0. 342) 

LMSC 

101 

73. 6 

f 


2. 92 

- 

2. 3 


0. 112 (0. 367) 

MDAC 

100 

75. 0 

b 


3. 0 


0 

150 x 150 

0. 0866 (0. 284) 

LMSC 

122 

100. 4 

f 


3. 3 

+ 

10. 0 

30 x 160 

0. 0750 (0. 246) 

LMSC 

142 

100 

b. 

f 

2. 82 

- 

6. 0 





120 

d 


3. 38 

+ 

12. 7 

24 x 110 

0. 0634 (0. 208) 

LMSC 

168 

138 

f 


3. 28 

+ 

9. 3 

*a. Kressilk 

b. Wintec c. 

Cambridge 

d. Jelliff e. Capital Westward f. LMSC 
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or 


“-•Prte 

The term in parentheses is a function of screen configuration, wire size, 
weave, etc. and is a direct indication of screen flow loss for laminar flow. 
Similarly, for turbulent flow, with R large, a / R is much smaller than |3, and 

f _ AP t 2 p §c _ q 
Q b pV2 


or 

AP = (5) 

Le DJ “c 

Again, the term in parentheses indicates the magnitude of the pressure loss 
for turbulent flow. 

The screen survey performed determined the important dimensional and 
geometric qualities of the screens, together with the flow-loss indicators 
specified by equations (4) and (5). The survey was limited to screen 
materials with application to aerospace needs, and thus started with 20 -mesh 
( ~ 1, 000 micron pore size) screens. (There are hundreds of meshes coarser 
than 20 mesh, made by dozens of suppliers, which are not included in the 
survey. ) 

It should be noted that screen nomenclature is based on English 
Engineering units, (e. g. mesh count in wires per inch, wire gage diameters 
in inches, etc. ) and no attempt was made to convert the screen survey or 
nomenclature to the Systeme International (SI) units because of the likelihood 
of confusion. Conversion factors for the reader's convenience are included in 
the results of the screen survey shown in Appendix A. The screens are 
organized as to weave: square, twilled square, plain Dutch, reverse plain 

Dutch, and twilled Dutch. Only woven screens are included; sintered or 
calendered mesh are not shown (see Appendix A for wire-cloth terminology). 
The principal suppliers of the full spectrum of screens are tabulated 
(although there are many other suppliers for some of the common screens) 
and coded in the survey. Most of the screens shown can be made from any 
metallic material which can be supplied as wire; all the screens shown are 
available in stainless steel. The weights and costs given are for stainless 
steel; for other materials, the weights and costs can be found from the 
supplier. Only representative (and approximate) costs are shown, since the 
costs do not increase appreciably until the very fine micronic sizes are 
encountered. For the Dutch weave screens, the manufacturer's rated 
absolute pore size in microns is shown (if known). The twilled Dutch screens 
shown on the third page of Appendix A are all the products of the Unique Wire 
Weaving Co. Inc. and are the only micronic grade of twilled Dutch screens 
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woven in the USA (all others are imported, generally from Germany, 
Switzerland, or Japan). The wire diameters used for weaving these screens 
are considered as proprietary information by Unique Wire Weaving, Inc. and, 
therefore, the geometric data and flow loss parameters are not shown, 
although the absolute screen ratings, based on alcohol bubble-point tests, are 
shown. It will be noted that the cost of these screens is two to four times 
that of similar imported screens. 

All of the screens shown in the survey, in stainless steel, are available 
from stock except the finest twilled Dutch meshes. All are also available, in 
virtually any material, on special order. The overall screen size, width, 
and length varies between meshes and manufacturers and should be checked 
on an individual basis with the supplier. 


Screen Selection 

Figure 4, shown previously, indicated the range of experimental LH 2 - 
tested pore sizes (bubble points) for screens, together with the largest tank 
dimension for the six tanks specified (see Table I). To facilitate the Task 3 
analysis, a realistic spectrum of pore diameters were chosen. 

Two observations are apparent from Figure 4: first, if 0. 1 g is the only 
selection criteria, then the circular pore diameter need not exceed the range 
of 10p to 50p (a span of less than an order-of-magnitude); second, based on 
previous LH 2 data, there is a data gap from 20pto 40 (j.. 

To obtain a pore size spectrum that will assure adequate coverage of the 
conditions of the Task 3 analysis, it was recommended that an acceleration 
range of 0. 1 g to 0. 0 1 g be used as the screen selection criterion - giving a 
pore size range of lOpto 500p. This range allowed selection of different 
weaves as well as coverage of the complete range of interest for the Task 3 
analysis . 

Selection of screens in this pore size range which will give the best 
potential performance makes it necessary to obtain screens with the maxi- 
mum ratio of bubble-point -to-flow-loss pressure drop for a given bubble 
point. Since the bubble point varies^as l/p a [see equation (2)] and the 

laminar flow loss varies as a m , the performance ratio can be com- 
puted for each screen, as shown plotted versus 1/ p a in Figure 6. This fig- 
ure, together with other selection criteria, such as cost, enabled the 
selection of the higher performance screens for use in this program. Based 
solely on Figure 6, the screens selected would have included the following: 
450 x 2,750, 325 x 2, 300, 850 x 850, 635 x 635, 500 x 500, 165 x 800, 

325 x 325, 50 x 250, 150 x 150, and 24 x 110 mesh. However, the 

450 x 2, 7 50, 850 x 850, and 63 5 x 635 screens have many practical problems 

associated with their use. These screens are very expensive (about 

$1, 400/m 2 ) and it is difficult to obtain large areas of screen with consistent 

performance, i. e. , without flaws or bubble point reduction over the entire 

sheet. The 850 x 850 and 635 x 635 mesh are woven of extremely fine wire 

and are very flimsy; it is difficult to visualize an actual WSL constructed of 
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Figure 6. Bubble-Point Flow Loss Ratio Parameter versus Bubble Point 



such flimsy material, and capable of tank installation in a reliable fashion. 
Finally, bubble-point data in LH;? were not available for the 850 x 850 and 
63 5 x 635 mesh, so that it was not certain that these screens would fill the 
gap in Figure 4 between 20 (jl and 40p. 

Therefore the 450 x 2, 750, 850 x 850, and 635 x 635 screens were not 
selected. Rather, the 200 x 1, 400 mesh and the 720 x 140 reverse Dutch 
were used together with the 500 x 500 mesh. This gave two distinctly dif- 
ferent weaves with about the same bubble point for Task 2 evaluation. The 
final screen selection to cover the entire range from 10pto 50 Ojjl with samples 
of every type of weave and in a cost-effective fashion is shown in Table III. 
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TABLE III. - FINAL SELECTION OF SCREENS 


Mesh 

Wire diameter, 
in. 

Weave 

Experimental 
circular 
pore size, 
microns 

Source 

Approximate 

cost 

$/m2 

325 x 2, 300 

0. 0015/0. 001 

Twilled Dutch 

12 

TET/Kressilk 

405 

200 x 1,400 

0. 0028/0. 0016 

Twilled Dutch 

19 

Gerard Daniel 

161 

500 x 500 

0. 001 

Twilled Square 

34 a 

TET /Kressilk 

389 

720 x 140 

0. 0014/0. 0043 

Reverse Dutch 

35. 1 

TET/Kressilk 

116 

165 x 800 

0. 0029/0. 002 

Twilled Dutch 

44. 5 

Gerard Daniel 

129 

50 x 250 

0. 0055/0. 0045 

Plain Dutch 

86 

Gerard Daniel 

55 

150 x 150 

0. 0026 

Square 

138 

Jelliff 

55 

24 x 110 

0, 015/0. 0105 

Plain Dutch 

168 

TET/Kressilk 

40 

60 x 60 

0. 0075 

Square 

310 a 

Jelliff 

35 

40 x 40 

0. 01 

Square 

510 a 

Jelliff 

31 

a Estimated 



















DETERMINATION OF SCREEN CHARACTERISTICS 


Bubble -Point Determination 

Screen samples were procured from the sources of the selected screens 
shown in Table III. The samples ranged from 0. 3 m to 0. 46 m by 1. 22 m 
(12 to 18 in. by 48 in. ) with the shute wires in the long direction. All exper- 
imental test specimens were fabricated from the single sample of each 
screen. The bubble-point specimens were 3 -cm circles and the flow-through 
specimens were 6. 5-cm circles cut from the same general area of the screen 
sample. 

The bubble-point and flow-through specimens only were cleaned, using 
cleaning procedures developed especially for this study after consultation 
with Wintec, Inc. , a leading firm in the fields of filter fabrication, cleaning 
of filter materials, and handling of fine-mesh screens. It was determined 
that the polyurethane adhesive which would be used to fasten the bubble -point 
specimen to the holding fixture was not compatible with long-term exposure 
to trichloroethylene, detergent, water, or isopropyl alcohol. MDAC experi- 
ence showed no problems with the adhesive during short term isopropyl 
alcohol bubble-point tests. However, it was decided to limit the screen 
cleaning material to Freon PCA with which the polyurethane adhesive is 
compatible. 

It was determined that the adhesive was unaffected by sonic cleaning in 
Freon PCA, and that there was no requirement for a particle count following 
cleaning for the bubble-point specimens. Since the flow-through specimens 
are oriented to the LH2 flow with the finest screen upstream, particle counts 
served no purpose because particles can go through the coarser downstream 
screens. The final NASA-approved screen cleaning procedure is shown in 
Table IV. The screen bubble-point specimens were precleaned in nitric acid 
(step 1 in Table IV), bonded to the specimen holders with EC3901 primer 
(3M Company) and Sta-Bond U-135 polyurethane adhesive per the bonding 
procedures of (Appendix A, ref. 11), then air-cured for 24 hours, and 
followed by an oven-cure at 71°C (160°F) for 24 hours. The specimens were 
bubble-point checked using ACS Reagent Grade isopropyl alcohol and helium, 
and then cleaned per the remaining steps of the procedure of Table IV. The 
bubble-point specimens (stainless steel elbows with the screen samples 
bonded to them) were arranged as shown in Figure 7 and suspended within a 
0. 17-m^ (45-gallon) LH 2 Dewar. The 10 samples were mounted so that they 
were visible through windows at the bottom of the Dewar. 

An LH2-level scale, indexed to the center of the bottom row of samples, 
was used to provide LH 2 head corrections to the individual bubble points, 
depending on where breakdown occurred. A 1389-G (ice point) platinum 
resistance temperature sensor was used to obtain the proper LH 2 temper- 
ature of 25. 2°K (45. 4°R), corresponding to saturation at 34. 5 N/cm^ 

(50 psia), in the vicinity of the samples. Dewar pressure control maintained 
the proper saturation conditions. Pressurization of the bubble-point speci- 
mens with gaseous hydrogen was accomplished through individual needle 
valves in the pressurization manifold. 
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TABLE IV. - SCREEN CLEANING PROCEDURE 


NOTE: All liquids and gases to be filtered through 0.45 micron membrane 
filter. 


1. Soak test item in pure nitric acid for 10 minutes minimum. 

2. Soak test item in Freon PC A for 15 minutes minimum. 

3. Drain Freon and immerse test item in a beaker of fresh Freon PCA 
and expose to an ultrasonic field of 16 to 20 Hz and a minimum 
intensity of 0.5 watts per square cm of tank bottom of a period of 

3 to 4 minutes. 

4. Rinse with Freon PCA for 1 /2 to 1 minute. 

5. Repeat steps 3 and 4. 

6. Purge dry with 0.45 micron filtered GN 2 . 

7. Package in clean Nylon 6 film. 

8. Overbag in 6 mil polyethylene. 


Installation of the Dewar into the LH 2 facility is shown in Figure 8. The 
Dewar was filled to about 0. 1 m (4 in. ) above the screen specimens (see 
Figure 7) and allowed to saturate at 34. 5 N/cm^ (50 psia). The LH 2 temper- 
ature in the vicinity of the screen samples was controlled to 25. 2 ± 0. 03°K 
(45. 4 ± 0. 05 °R) by controlling the Dewar pressure. The LH 2 was allowed to 
boil down to about 1/2 cm above the top row of screen specimens to minimize 
head effects on bubble point. The appropriate valve in the pressurization 
manifold was opened, and GH2 was slowly bled to the sample through the GH 2 
metering valve. The gas pressure behind the sample was slowly increased 
and monitored by a Merriam 150 -mm manometer and a precision hook gauge 
(Microtector) until bubbles appeared, as viewed through the viewports. The 
LH2 level and the point of bubble emergence on the sample were recorded. 

The primary bubble-point pressure measurements were made with the 
Microtector, which after several trials, gave repeatable results within less 
than 0. 025-mm H 2 O. The 150-mm manometer reading was also recorded, 
as a backup, and agreed with the Microtector value within less than 0. 5-mm 
H 2 O. Following the LH 2 bubble-point tests, the screen samples were again 
bubble-point checked with ACS Reagent Grade isopropyl alcohol. The surface 
tension of the alcohol used for the bubble-point tests was measured with a 
DeNuoy tensiometer. The average of several trials gave a value of 23.8 
dynes /cm at 22.2°C (72 °F), which is somewhat in excess of the standard 
value. It is possible that water absorption by the rather hygroscopic 
alcohol caused the surface tension increase. The alcohol and LH 2 bubble - 
point data (corrected for head effect) are summarized in Table V. The 
alcohol data are shown corrected to the standard isopropanol surface tension 
value of 21. 15 dynes/cm at 25°C (77°F), and also expressed in meters of 
34.5 N/cm.2 (feet of 50 psia) LH^, using the properties given in Reference 6. 
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Figure 8. LH2 Bubble Point Test Setup 



TABLE V. - ALCOHOL AND 34. 5 N/cm (50 psia) LH BUBBLE-POINT TEST DATA 


Screen 
(wire dia) 


Pre LH 2 test alcohol data 


Post LKL, test alcohol data 


Measured Predicted bubble Measured Predicted bubble 
bubble point point in meters of bubble point point in meters of 
in inches H O 34. 5 N/cm 2 (feet in inches H ? 0 34 - 5 N/cm z (feet 
(standard) of 50 psia) LH 2 (standard) of 50 psia) 


LH 2 data 

Bubble point 
in meters of 
34. 5 N/cm 2 (feet 
of 50 psia) LH 2 


325 x 2, 300 
(0.0015/0.001) 

23. 90 

0. 5090 

(1. 670) 

200 x 1,400 

(0. 0028/0. 0016) 

16. 55 

0. 3523 
(1.156) 

500 x 500 
(0. 001) 

7. 73 

0. 1646 
(0.540) 

720 x 140 

(0. 0014/0.0043) 

7. 82 

0. 1664 
(0. 546) 

165 x 800 
(0. 0029/0. 0020) 

6. 70 

0. 1426 
(0.468) 

50 x 250 
(0.0055/0.0045) 

3.45 

0. 07346 
(0. 241) 

150 x 150 
(0. 0026) 

2. 25 

0. 0480 
(0. 1575) 

24 x 110 

(0. 015/0. 0105) 

1.685 

0. 0359 
(0. 1178) 

60 x 60 
(0. 0075) 

0. 998 

0. 02121 
(0. 0696) 

40 x 40 
(0.010) 

0. 620 

0. 01323 
(0. 0434) 




0.4602 
( 1.510) 

0.4815 

(1.580) 

0.3517 

(1.154) 

0. 3377 
(1.108) 

0.1652 
(0. 542) 

(0. 1704) 
0. 146t> 

/ (0. 559) a 
\ 0. 479 b - 

0.1664 

(0.546) 

0. 1767 
(0. 580) 

0. 1436 
(0.471) 

0. 1228 
(0.403) 

0. 0719 

0. 0682 


0. 04816 

0. 0460 

(0. 1580) 

(0. 1510) 

0.03612 

0. 0336 

(0.1185) 

(0. 1105) 

0.02118 

0. 0230 

(0. 0695) 

(0. 0754) 

0. 01332 

0. 0170 

(0.0437) 

(0.0559) 


Manufacturer's 
Rated Absolute 
pore size 
(microns) 




General breakdown 
Single wire breakdown 







The actual LH 2 bubble-point data are compared to the predicted value 
based on the alcohol data in Figure 9, and to the absolute screen pore size 
rating in Figure 10. From Figure 9, the only sample that showed a differ- 
ence in the pre-LH^ test and post-test alcohol data was the 325 x 2, 300 
screen, and this unexplained difference was less than 10%. 

The 500 x 500 twilled square screen exhibited anomalous behavior in the 
LH 2 bubble-point tests. Bubble retention failure occurred at a lower value 
than expected, and was confined to a single horizontal wire location, which 
would be easily seen with the naked eye. This single wire location was not 
the failure point in either the pre- or post-test alcohol bubble-point check, 
where failure occurred about 1 cm above the horizontal wire location (which 
could still be seen, but was not a failure point until general screen breakdown 
was induced). During the LH 2 tests, general screen breakdown occurred at 
a value closer to that predicted (see Table Y). The reason for this anomaly 
may be that the 500 x 500 screen is extremely flimsy, and cryogenic con- 
traction may have somehow caused the wire location to open. 

The 40 x 40 square weave screen also showed anomalous behavior in the 
LH 2 tests by exhibiting a much higher bubble point than predicted from the 
alcohol data. As can be seen from Table V, this screen has a very small 
bubble -point value in LH 2 , and thus the effects of all sources of error (such 
as instrument zero shift, head contribution, etc. ) are magnified, even though 
great care was taken to acquire accurate data. 

Many careful repetitions of the bubble -point measurements, with 
instrument zero rechecks between them, were made with the 40 x 40 screen, 
with the same unexplained result. It was thought that the vertical orientation 
of the screen might have influenced the bubble-point value, so during the 
post-test alcohol checks, the 40 x 40 screen and the 60 x 60 screen were also 
alcohol tested in a horizontal orientation facing upwards. The bubble-points 
in both orientations were virtually identical, and in fact the horizontal orien- 
tation had a higher bubble -point by about 5% for both screens. (This can 
perhaps be explained by gas interface distortion due to buoyancy just before 
bubble -through which would tend to decrease the bubble point in the vertical 
orientation. ) 

The superior bubble-point performance of the 40 x 40 screen is also 
evident in Figure 10, compared to the other screens. The line in Figure 10 
is for a <j>' = 2. 6 in the bubble-point equation (1). 

h 

gP D 

where D is the manufacturer's rated absolute pore size (inscribed circle 
of the largest pore in the screen). The value of 4>' is 4. 0 for circular pores, 
however, in general, the value of cf>' is degraded for Dutch twill and square 
weave screens which have triangular and square pores, respectively. In 
previous bubble -point work with saturated LH 2 at 10. 1 N/cm^ (14. 7 psia) 
Table II), <j>' was found to be more nearly equal to 3. 0, while for our tests the 
value of (j) 1 is 13% lower. It is perhaps coincidental that the heat of vaporiza- 
tion at 34. 5 N/ cm^ (50 psia) is 10% lower than that at 10. 1 N/cm^ (14. 7 psia). 
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but heat transfer to saturated liquids may have a noticeable effect on bubble 
point. It was noticed during bubble -point testing that lights, directed through 
the viewports onto the screens, caused fine bubbles to emanate from the 
screens, especially the fine meshes. Therefore, all bubble-point observa- 
tions were made with rather dim available light which fortunately was ade- 
quate for viewing. From Figure 10, all of the screens follow the <J>' line 
rather well, except for the three coarsest screens, which gave higher per- 
formance and 4»* values closer to 3. 0 (see Table V). It may be significant 
that the wires which make up the three coarsest screens have cross-sectional 
(heat conduction) areas at least nine times that of the finer screens. This, 
combined with fewer wire joint or pore nucleation sources, may contribute 
to reducing the effects of heat transfer, if any, resulting in increased per- 
formance. In any event, it would be accurate or conservative to use the 4>' 
value of 2. 6 in the bubble-point equation for a functional description of 
bubble-point, based on the manufacturer's absolute rating for screens in 
34. 5 N/cm2 (50 psia) saturated LH 2 . 


Flow Test Apparatus Design and Fabrication 

To characterize the screen flow-through loss and channel flow loss 
parameters, all 10 screens had to be tested at several flow-through rates 
and, for the channel flow loss tests, at several flowrates and channel 
spacings. This required obtaining literally hundreds of data points and 
design of a test apparatus to efficiently determine these data accurately but 
with minimum apparatus modifications represented a significant challenge. 

A schematic representation of the apparatus is shown in Figure 11. The 
basic channel flow apparatus was a box formed of a number of channels in 
series, as shown. A movable partition lined with screen formed one side of 
two channels in series, so that one partition position gave data on two annulus 
spacings. The LH 2 flow entered through a long entrance tube and manifold to 
minimize entrance effects. The flow-through specimens were integrated into 
this inlet line and provided flow-through loss data simultaneously with the 
channel loss data. The entire apparatus was installed in a large insulated 
pipe which served as a pressure vessel for the 34. 5 N/cm^ (50 psia) LH 2 . 

The LH 2 , after leaving the channel box, flooded the interior of the pipe and 
provided thermal isolation before leaving the pipe through the back pressure 
system. The initial design problem was to place the screens in channels with 
gaps that were representative of the tank application of these screens, so 
that the data obtained would be usable in the Task 3 analysis without excessive 
extrapolation. The tankage systems specified are shown in Table VI, with 
the annulus gap for each tank for an annular residual of 1% to 5% of tank vol- 
ume shown. In all cases, the apparatus gap range was chosen to span 
the annulus gap range shown. From the static retention analysis 
described previously for the six tanks over a range of 0. 1 to 0. 01 g's, 
the required screen bubble point in microns for static retention of the 
head imposed by the largest tank dimension, L, for 34.5 N/cm^ 

(50 psia) LH 2 is shown. Finally, the screens selected for each channel pass 
are shown in the last column. The finest screen, 325 x 2, 300, was used with 
the finest static retention requirement. For direct AP comparison between 
twilled Dutch and twilled square weaves, the 500 x 500, which has the same 
shute wire diameter (assumed to be an independent variable for pressure 
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Figure 1 1 . Flow-Loss Test Apparatus 




TABLE VI. 



Vol 

m 3 (ft 3 ) 

L»/D 

L 

m(ft) 

Annulus 

l%-5% 

gap, 

cm (in. ) 

1 

141. 6 
(5, 000) 

4 

14. 65 
(48) 

0. 84-4. 19 
(0. 33-1. 65) 

2 

14. 16 
(500) 

4 

6. 80 
(22.3) 

0. 38-1. 96 
(0. 15-0. 77) 

3 

14. 16 
(500) 

2 

4. 42 
(14. 5) 

0. 46-2. 31 
(0. 18-0. 91) 

4 

14. 16 
(500) 

1 

3. 00 
(9.9) 

0. 51-2. 52 
(0. 20-0. 99) 

5 

1. 416 
(50) 

2 

2. 06 
(6.7) 

0. 213-1. 069 
(0. 084-0. 421) 


0. 231-1. 161 
(0. 091-0. 457) 






loss), was also used in the same channel. The next finest screen, 

200 x 1,400, was in the channel with the next finest static retention require- 
ment; for direct pressure loss comparison, a reverse Dutch weave was 
used. The coarsest mesh, 40 x 40, was used in the channel with the 
coarsest static retention requirement. A plain Dutch weave (24 x 110) of 
nearly the same wire diameter was used in the same channel for the pressure 
loss comparison. The remaining four screens were placed in the remaining 
channels in generally increasing micron size. The configurations 1 and 2 
noted in Table VI are shown schematically in Figures 12 and 13, with 
only the vertical dimensions approximately to scale. Configuration 1 had 
six passes and six screens mounted as shown in the channels, plus the six 
Dutch weave screen flow-through specimens. Configuration 2 had fixed 
partition A (see Figure 12) removed to give the large gaps necessary for the 
325 x 2,300 and 500 x 500 screens. The remaining two screens were mounted 
in the channels of Configuration 2, along with the four square weave flow- 
through specimens. The channel flow test matrix is shown in Table VII. The 
as -built channel spacing and planned volumetric flowrates are shown. The 
actual flowrate obtained during testing varied continuously but approximated 
the values given in Table VII. 

The apparatus was made with all-welded construction. The critical 
design problem was to configure the sides of the apparatus to provide the 
close tolerances required on the channel height over the 2.26-m (89-in.) - 
long apparatus. Machining the grooves was rejected because of problems of 
extreme length and alignment, and requirements for thick material which 
leads to potential welding problems and warpage. Instead, the grooves were 
accurately rolled into the thin 0.5 -mm (0. 020-in.) material. This minimized 
the heating necessary for welding and minimized weld damage and warpage. 
Some difficulty in rolling caused by work-hardening resulted in the as-built 
channel gaps being somewhat different from the planned gaps (Tables VI and 
VII). A large number of static pressure sensing taps were used to determine 
the pressure loss along the channel. The pressure taps were arranged to 
provide data in both screen sample partition positions. The length between 
the pressure taps was arranged to give three length values (e.g. , 1 to 2, 2 
to 3, and 1 to 3), which allowed a check on the linearity of pressure loss 
with length. The pressure taps were situated alternately on opposite sides of 
the apparatus to avoid interference and allow full use of the pressure-vessel 
flange for external routing of pressure sensing lines. The pressure taps were 
also arranged to give a channel entrance L/D ratio of at least 5 to allow flow 
smoothing leaving the entrance manifold and the 180 degree "mitered" bends. 
The entrance manifold contained vanes to equally distribute the flow from 
the entrance line to the apparatus. The assembled flow loss test apparatus is 
shown in Figure 14. The pressure sensing lines were routed horizontally 
(to avoid head errors), passed through the pressure vessel flange, and 
terminated in the needle-valve pressure-tap panel. The principal pressure 
instruments used were a pair of Dwyer "microtector " electronic hook gauges 
with a range of 0 to 5.08 cm (0 to 2 in.) of H^O and a sensitivity of 0.001 cm 
(0. 0005 in. ) of H^O. A Merriam 0 to 7 6 -cm (0 to 30 in. ) of H^O manometer 
was also used for the higher pressure loss measurements. A complete 
pressure loss analysis through the entire flow apparatus was performed 
to determine the necessary flowhead and ensure that at the nominal flow- 
rates (Table VII) the pressure drops would be within the range of the 
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TABLE VII. 


CHANNEL FLOW TEST MATRIX 


Channel Width = 12. 5 cm (4. 94 in. ) 


Screen 

Actual channel 
spacing depth 

Approximate 
equivalent 
annulus gap, 
% 

Channel 
length, 
cm (in. ) 

Flow Rate, 
mV sec (gpm) 

cm 

in. 

40 x 40 

0. 455 

0. 179 

1. 95 

33 66 99 

0. 001, 0. 0015, 0. 002 


0. 554 

0. 218 

2. 4 

(13 26 39) 

(16, 

24, 32) 


1. 034 

0. 407 

4. 5 






1. 133 

0. 446 

4. 9 





24 x 110 

0. 457 

0. 180 

2. 0 






0. 531 

0. 209 

2. 3 






1. Oil 

0. 398 

4. 4 






1. 085 

0. 427 

4. 7 





720 x 140 

0. 622 

0. 245 

1. 6 






0. 747 

0. 294 

1. 95 






1. 514 

0. 596 

4. 0 






1. 638 

0. 645 

4. 3 





200 x 1, 400 

0. 635 

0. 250 

1. 7 






0. 762 

0. 300 

2. 0 






1.499 

0. 590 

3. 9 






1. 626 

0. 640 

4. 3 





60 x 60 

0. 587 

0. 231 

1. 15 






1. 252 

0. 493 

2. 5 






1. 948 

0. 767 

3. 8 






2. 614 

1. 029 

5. 1 





50 x 250 

0. 589 

0. 232 

1. 15 






1. 234 

0. 486 

2. 4 






1. 963 

0. 773 

3. 9 



0. 001, 0. 0015, 0. 002 


2. 609 

1. 027 

5. 1 



(16, 

24, 32) 

150 x 150 

0.612 

0. 241 

1. 35 



0.003, 0.0035, 0.004 


1. 278 

0. 503 

2. 8 



(48, 

56, 64) 


1. 974 

0. 777 

4. 3 






2. 639 

1. 039 

5. 8 





165 x 800 

0. 610 

0. 240 

1. 35 






1. 255 

0. 494 

2. 7 






1. 984 

0. 781 

4. 3 

33 66 99 




3. 629 

1. 035 

5. 75 

(13 26 39) 



500 x 500 

1. 090 

0. 429 

1. 3 

99 132 198) 




1. 643 

0. 647 

1. 95 

(39 52 78) 




(2. 890) a 

(1. 138) a 

3. 45 






(3. 444) 

(1. 356) 

4. 1 





325 x 2, 300 

(1. 090) 

(0. 429) 

1. 3 






(1. 643) 

(0. 647) 

1.95 






2. 890 

1. 138 

3. 45 

99 132 198 

0. 003, 0. 0035, 0. 004 


3. 444 

1. 356 

4. 1 

(39 52 78) 

(48, 

56, 64) 


Parentheses indicate optional measurement if required 
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Figure 12. Configuration One Flow Apparatus (Vertical Dimensions Approximately to Scale) 













instrumentation, and more important, was in the appropriate range for later 
use of the data in the Task 3 analysis. The flow-loss test apparatus was 
inserted in the pressure vessel and installed in the MDAC LH 2 test facility 
as shown in Figure 15. The LH 2 was allowed to saturate in the storage tank 
at 34. 5 N/cm.2 (50 psia) prior to use. During testing, LH 2 at 34. 5 N/cm^ 

(50 psia) flowed through a 10 -micron filter, through the inlet control valve, 
the flowmeter, and the specimens in the apparatus. The flow then emptied 
into the pressure vessel where it provided thermal shielding of the specimens 
and was exhausted through another control valve into a 3. 79-m^ ( 1000 -gallon) 
catch tank. When the catch tank was full, testing was halted and the LH 2 was 
returned to the storage tank for reuse. Two previously anticipated problems 
were encountered: there was some difficulty in obtaining single phase LH 2 

flow through the specimens. This was overcome by increasing the storage 
tank pressure by up to 3. 5 N/cm^ (5 psi). In addition, it was found that a 
minimum of about 0. 001 mV sec (16 gpm) of LH 2 flow was required to ensure 
single-phase flow through the specimens. The presence of two-phase flow 
was detected by low pressure drop across the screens, and by cyclic surging 
of the screen pressure drop. The flowrate was determined using a turbine 
type flowmeter, and the LH 2 temperature entering and leaving the apparatus 
was determined using platinum resistance sensors. These data, together 
with time, were recorded on a digital data paper tape system. The time was 
also manually recorded together with the pressure-drop data for later 
correlation with the paper-tape data. 


Screen Flow-Through Test Results 

The screen flow-through specimens were welded between stainles s -steel 
washers, and sealed between spacers in the LH 2 inlet line to the flow test 
apparatus. The screen area exposed to the flow was 23 cm2 (3. 56 in. 2) anc L 
was the same as the inside area of the flow passage. The entrance section in 
the inlet line up to the first specimen was 1. 4 m (55 in. ) long (L/D = 26) 
which is believed adequate for inlet flow smoothing. The pressure drop 
across the screen, H, was determined in inches of H 2 O and converted to 
meters of saturated LH 2 at 34. 5 N/cm.2 (50 psia). The flowrate together with 
the screen (flow passage) area determined the fluid approach velocity, Y, in 
m/sec. As described in the section on Screen Survey, the pressure drop for 
flow through a screen can be described in terms of a friction factor, f, and a 
Reynold's number, R, in the manner of Armour and Cannon (ref. 10). The 
correlation is 


f = 



( 6 ) 


where 


H « 2 Dg c 
f " V 2 Qb 


and 


R 


P 


pa 
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Figure 15. LH 2 Flow Loss Test Setup 



(See Symbols) and a and P are experimentally determined constants. 

Substitution of f and R gives the following expression relating pressure drop, 
H, and fluid approach velocity, V: 


H = P 


" ^ 2 


’ Qb 

Qba 

pV + p 

2 p 

£ Sc 

e 2 Dg 

to c 


V' 


H = AV + BY 


(7) 


From equations (6) and (7), it can be seen that in the laminar flow regime 
(small V and R), f is essentially only a function of R, and H depends essen- 
tially only on V, while in the turbulent flow regime, f is essentially constant, 
and H essentially depends only on V^. For the LH^ tests, a minimum flow- 
rate of about 0.001 m^/sec (16 gpm) was required to maintain good quality 
LH 2 in the apparatus. This, together with the low viscosity and (relatively) 
high density of the LH 2 , resulted in the LH 2 test data being in the turbulent 
Reynold's number regime only. To evaluate the functional constants in the 
laminar regime, additional tests were made with both the channel apparatus 
and with another flow apparatus used for the Reference 12 flow-loss tests, 
using ambient gaseous nitrogen as the flow medium. The relevant screen 
characteristics and the experimentally determined a , p, A and B for the 10 
selected screens are shown in Table VIII. The data, plotted in Armour and 
Cannon form, are shown in Appendix B for all screens, with typical 
examples shown in Figures 16 and 17. 


Figure 16 shows the available data for the 325 x 2, 300 screen, including 
gas data from the MDAC-MSFC contract NAS8-27685 (ref. 12) which agree 
with our data rather well. Water test data from Kressilk (ref. 13) and 
Wintec Corp. (ref. 14) are also shown, together with the generalized Armour 
and Cannon correlation for all screens, which overpredicts the friction 
factor and thus the pressure loss for this screen by 170 percent. 


Also shown in Appendix B are data from all other known sources, 
i. e. , GDA (ref. 15), MDAC IRAD (ref. 16), NAR (ref. 17), and 
MDAC-MSFC contract NAS8-27685 (ref. 18). These data have been 
normalized to agree with the geometric characteristics of our screens 
(especially specified pore diameter) as shown in Table VIII. 


The correlations for the square weave screens are all shown in Fig- 
ure 17, which exhibits a number of interesting aspects. First, in the laminar 
regime, all four screens are reasonably well represented by a single corre- 
lation value of a. Second, in the turbulent regime, assuming one -dimensional 
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TABLE VIII. - FLOW LOSS CORRELATION 


o 


Screen 

D, 

m (ft) 

a, 

l/m ( 1/ft) 

b, 

m (ft) 

€ 

Armour and 
Cannon Form 

f -t + (3 

H = AV + BY 2 

H in m of 
34. 5 N/cm 2 LH 2 
V in m/ sec 

H in ft of 
50 psia LH2 
V in ft/ sec 

a 

P 

A 

B 

A 

B 

325 x 2, 300 

0. 000005 

110, 235 

0. 000089 

0. 245 

3. 2 

0. 19 

1. 14 

0. 6919 

1. 14 

2. 27 


(0.0000164) 

(33, 598) 

(0. 000292) 








200 x 1, 400 

0. 000010 

65, 390 

0. 0001524 

0. 248 

4. 2 

0. 20 

0. 885 

0. 6126 

0. 885 

2. 01 


(0. 0000328) 

(19, 930) 

(0. 0005) 








720 x 140 

0. 000015 

32, 954 

0. 0001804 

0. 514 

11. 0 

0. 47 

0. 162 

0. 2627 

0. 162 

0. 862 


(0. 0000492) 

(10, 044) 

(0. 000592) 








165 x 800 

0. 000025 

41, 360 

0. 0001753 

0.426 

3. 3 

0. 17 

0. 108 

0. 0805 

0. 108 

0. 264 


(0. 000082) 

(12, 606) 

(0. 000575) 








50 x 250 

0. 00005 

13, 075 

0. 000369 

0. 611 

13. 5 

0. 26 

0. 045 

0. 0631 

0. 045 

0. 207 


(0.000164) 

(3, 985) 

(0. 00121) 








24 x 110 

0. 000115 

5, 889 

0. 000914 

0. 572 

8. 61 

0. 52 

0. 0165 

0. 1554 

0. 0165 

0. 51 


(0. 000377) 

(1, 795) 

(0. 003) 








500 x 500 

0. 0000254 

65, 495 

0. 0000509 

0. 584 

5. 7 

0. 65 (0. 77) a 

0. 0554 

0. 03298 

0. 0554 

0. 1082 


(0.0000833) 

(19,962) 

(0. 000167) 





(0. 04267) a 


(0. 140) a 

150 x 150 

0. 000103 

19, 916 

0. 000132 

0. 671 

5. 7 

0. 50 (0. 50) a 

0. 0101 

0. 01347 

0. 0101 

0. 0442 


(0. 000339) 

( (6, 070) 

(0. 000433) 





(0. 01347) a 


(0. 0442) a 

60 x 60 

0.000233 

(8, 137) 

0. 000381 

0. 612 

5. 7 

0. 40 (0. 61) a 

0. 00585 

0. 01820 

0. 00585 

0. 0597 


(0. 000764) 

(2, 480) 

(0. 00125) 





(0. 02761) a 


(0. 0906) a 

40 x 40 

0.000381 

5, 328 

0. 000509 

0. 662 

5.7 

0. 60 (0. 5 2) a 

0. 00287 

0. 01728 

0. 00287 

0. 0567 


(0. 00125) 

(1, 624) 

(0. 00167) 





(0. 01497) a 


(0. 0491) a 

Armour and 

Cannon (Refe 

rence 10) 



8. 61 

0. 52 






a Based on Ee = ^ ^ ^ 







incompressible flow with a sudden expansion, a potential flow analysis which 
balances momentum loss with pressure recovery behind the screen (ref. 19) 
results in the expression: 



( 8 ) 


where the Euler No. Eu = AP 2g c /pv2, and S is the screen solidity (fraction 
of closed area). This is equivalent to a constant friction factor where 


f 


Eu i^D 
2 Qb 



(9) 


which is a function only of screen geometric characteristics. 

The value of (3 and B based on equation (9) for our four square-weave 
screens is shown in Table VIII. There is reasonably good agreement for 
all of the screens, with the 40 x 40 screen value being perhaps a little high, 
and the 60 x 60 screen value being a little low. Shown for comparison is 
Wintec Corp. water data which shows very poor agreement, especially 
in the turbulent regime. The Wintec equation for the 60 x 60 screen has 
only the turbulent component which gives an equivalent f of 1. 82. 

Evaluation of the equivalent f for the ten 60 x 60 screens included in our 
screen survey (Appendix A) indicates that the equivalent f should be in 
the range of 0, 275 to 1. 07. This illustrates the potential unreliability 
of much available screen pressure loss data. 

2 

The data for pressure loss, in meters of 34. 5 N/cm LH 2 
( P = 64. 08 kg/m^) versus approach velocity in meters/sec are shown in 
Figure 18. The lines shown are equation (7), based on the a and |3 of 
equation (6) and shown in Table VIII for each of the 10 screens. The 
deviation of some of the data for the correlation indicates the desirability of 
obtaining the correlation over the entire laminar-turbulent flow regime 
based on dimensionless parameters such as f and R, rather than simply 
on pressure drop versus velocity data. It is believed that the correlations 
described in Table VIII adequately describe the screen flow-loss 
characteristics over the entire flow spectrum. 


Channel Flow Test Results 

The channel-flow screen specimens were screen samples 0. 125 m by 
1.0 m (5 in. by 40 in. ) (with the shute wires in the long direction of flow) 
bonded to stainless steel backup plates with the same polyurethane adhesive 
used for bonding of the bubble-point specimens. The adhesive was used 
sparingly and wiped from the top surface of the screens to ensure that the 
screens would exhibit the necessary roughness due to the weave. The com- 
plete matrix of test conditions shown in Table VII was performed. 


43 



PRESSURE DROP - FT OF 50 PSIA LH 2 


CR140 



Figure 18. Flow Loss Correlation - Pressure Drop Versus Velocity 
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PRESSURE DROP 


The data, plotted as head Loss in meters of 34. 5 N/crrm LH 2 versus chan- 
nel flow velocity in meters /sec with length (L) over channel height (s) as a 
parameter, are shown in Appendix B for all 10 screens, with a typical 
example shown in Figure 19. As anticipated in test planning, the pressure 
drop was linear with length: in Figure 19, the plain symbols are for the 
33 cm (13-inch) length, the primed symbols are for 66 cm (26-inch), and 
the double-primed symbols for the 99 cm (39-inch) length. Because the 
channel flow lies in the turbulent regime, the pressure drop varies 
essentially as the channel flow velocity squared, as shown by the slope of 
the lines on Figure 19. In general, the pressure drop increases with 
increasing L/s, but extreme channel height values cause jumps in the L/s 
correlation because of the effect of channel height on friction factor (Fig. 19). 


To evaluate this effect and define the important physical factors influ- 
encing the channel flow, a dimensionless analysis of the channel flow, in the 
manner of Moody (ref. 20) was performed. In pipe flow, the friction factor 
is defined by the Darcy formula (see Symbols): 


H f . H f D h w Z s 2 2g c 

L V 2 L Q 2 

D h 2 §c 


( 10 ) 


and at ordinary velocities f is a function only of Reynolds number, R, and a 
roughness parameter, e/D^, where 


P VD h _ 

p v 


(ID 


and e is a linear dimension representative of the absolute roughness of the 
surface. For non-circular pipes, in e/Dfr and equations (10) and (11) is 
the hydraulic diameter defined as four times the flow area divided by the 
wetted perimeter, or, for a rectangular channel of width w and height s. 


ws 

D, = 2 — ; — 
h w + s 


( 12 ) 


In Reference 20, Moody presents a graph giving the relationship 
between f, R, and e/Db- The pipe flow is characterized by four flow 
regimes: in the laminar flow regime, the flow loss correlation is of the form 


f 



(13) 


and the roughness has no effect on the flow, which is dominated by viscous 
effects. At a critical R of 2, 000 to 4, 000, ordinary laminar pipe flow makes 
a rather abrupt transition to turbulent flow. In this critical region, the f is 
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not well defined and this region is shown on the Moody graph as a rather wide 
band. Following the critical region, the flow enters the transition region, 
where f is a complex function of R and e/Dh. At high values of R, the flow 
enters the region of complete turbulence for rough pipes and R has no effect, 
so that 


f - f(e/D h ) (14) 

To determine the effects of e/D^ on our channel flow loss and the 
appropriate roughness parameter, our data were plotted on the Moody graph, 
as shown in Appendix B for all 10 screens, with typical examples shown in 
Figures 20 and 21. In these figures, the collection of data points at 
R = 105 or higher were the 34. 5 N/cm^ (50 psia) L,Hz data, while the 
single points at R = 5, 000 were data taken with ambient GNg. In our 
experiments, the direction of fluid flow was always in the direction of the 
shute wires, which should minimize the pressure drop due to the construction 
of the Dutch-weave screens. It was noted from the data that the f for the 
Dutch-weave screens was about half that for square-weave screens of about 
the same wire diameter and essentially identical geometry and flow 
conditions . 


To determine the value of e more precisely, it was computed, 
based on the experimental data, from the Colebrook equation (ref. 20) 
which describes the relation between f, R and e/Dh i- n the turbulent- 
transition regime: 



-2 log 



+ 


2. 51 \ 
R \/f / 


(15) 


The value of e based on the experimental data and actual channel spacing, 
and equation (15) is shown in Table IX, compared to e based on physical wire 
dimensions— i. e. , the wire diameter for the square-weave screens, and half 
the shute-wire diameter for the Dutch-weave screens. The values in 
Table IX substantiate the use of wire diameter (or half the shute-wire 
diameter) as the appropriate roughness dimension. It is, of course, desir- 
able that this roughness parameter be based on a general physical character- 
istic of the screen, rather than on specific experimental data for each 
screen. 

The physical rationale behind this is shown in Figure 22, which sche- 
matically shows a cross-section through the screens in the direction of flow. 
The construction of the Dutch-weave screens is such that the shute wires are 
pressed tightly together so that essentially only half the shute wire con- 
tributes "roughness." On the other hand, the square-weave screens are 
woven with definite "pitch", or spaces between wires, so that the flow sees 
the full wire diameter as "roughness." A similar effect wap found by 
Edwards and Sheriff (ref. 21), who determined that at low pitches ( 1 to 4 
diameters) the friction factor increases nearly linearly with pitch. For our 
case, the pitch of the square-weave screens is 2 to 2. 5 wire diameters, 
while for the Dutch-weave screens the pitch is effectively 1 diameter. 

It can be seen from the roughness values of Table IX, and from Fig- 
ures 20 and 21, that there is considerable data scatter. There are a number 
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VALUES OF (VO’) FOR WATER AT 60' F (VELOCITY In|5t.X OIAKETlR iN INCHES) 


0 10 


I 1 t f VALUES OF fVO’) FOR ATMOSPHERIC AlR AT 60'E 
2 4 I 6 8 lOl 20 I 40 1 60 I 1 ' 

r- ■ v ’! TfifiVi! 1 -- 1" i‘i rni 

trr. :r: i. til ! i" n : 


600 eooooo 


4 00 G OO 80 0 KX>0 

! eooo I 

4000 | 6 ooo 1 10,000 ; 

,l.i. . !ii , . , 


LA'JWtAP ! CRfTKAl TRANsmqtfv L-jJ-U J J. h 11 .1 : M .!.].] 1 !.: i ! J_J 11 

i FL0»V — +-ZD-N E — tr ZONE A A* •--4-4-!-+ COMPLETE TURBULENCE, ■ ROUGH ; PIPES- ■ - t - : ■■ — -~ 

\-- ttti nTrrnfr=-v : mt — r fri i ! ;ri m; ,-rrr: n 

J r ! 1" f rm^^Titi H rrff “ 4 “ r r ! r I 


£000 

•* 4000 fOpOSp OOl 

0-0,000 j I 

40,000 60,000 100,000 

.4* .<. r .A£...L. t -- 


linjltt 

' ~ lii 4 ri' It 

fri rinllt 

_ — ; ~rrm r — rr 

rrii ; nfrff 


rBsi 


150 X 150 


l0 3 2(I0 5 ) 3 4 5 6 8 IQ4 2(l0 4 ) 3 4 5 6 8 |q5 ZOO 3 ) 3 4 5 6 8 |q6 2U0 4 ) 3 4 5 6 8 |q7 > *^0')'*-3 _^5 6 8 |q8 

. - _ . .. _ \/r> / trr ctt2 \ 0 m -0Or> _ 7T**.£ 


REYNOLDS NUMBER R=-^_ (Vwg , D IN FT, V IN^g?) 


Figure 21. Dimensionless Channel Flow Loss Correlation — 150 x 150 
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TABLE IX. - SCREEN WIRE SIZE AND ROUGHNESS PARAMETERS 


U1 

o 


Screen 

Wire diameter - in. 
warp/ shute 

Roughness parameter - e - cm 
Based on experimental data 
and Colebrook equation 3 - 

Roughneso parameter - e - in. 
Based on experimental data 
and Colebrook equation 3 - 

Roughness parameter 
e - cm - (in. ) 
based on physical 
wire dimensions 

325 x 2, 300 

0. 0015/0.001 

0. 00127 
(2. 89)4 

0. 001295 
(3. 444) 



0. 00050 
(1. 1 38) c 

0. 00051 
(1. 356) 



0. 00127 (0. 0005) 

200 x 1, 400 

0. 0028/0. 0016 

0. 00353 
(0. 635) 

0. 00183 
(0. 762) 

0. 00239 
(1. 499) 

0. 00173 
(1. 626) 

0. 00139 
(0. 250) 

0. 00072 
(0. 300) 

0. 00094 
(0. 590) 

0. 00068 
(0. 640) 

0. 00203 (0. 0008) 

720 x 140 

0. 0014/0. 0043 

0. 00719 
(0. 622) 

0. 00462 
(0.747) 

0. 00617 
(1. 514) 

0. 00538 
(1. 638) 

0. 00283 
(0. 245) 

0. 00182 
(0. 294) 

0. 00243 
(0. 596) 

0. 00212 
(0. 645) 

0. 00546 (0.00215) 

165 x 800 

0. 0029/0. 002 

0. 00244 
(0. 610) 

0. 00312 
(1. 255) 

0. 00310 
(1. 984) 

0. 00241 
(2. 629) 

0. 00096 
(0. 240) 

0. 00123 
(0. 494) 

0. 00122 
(0. 781) 

0. 00095 
(1. 035) 

0. 00254 (0. 001) 

50 x 250 

0.0055/0.0045 

0. 00686 
(0. 589) 

0. 00584 
(1. 234) 

0. 00630 
(1. 963) 

0. 00460 
(2. 609) 

0. 0027 
(0. 232) 

0. 0023 
(0. 486) 

0. 00248 
(0. 773) 

0. 00181 
(1. 027) 

0. 00572 (0. 00225) 

24 x 110 

0.015/0. 0105 

0. 02159 
(0. 457) 

0. 01854 
(0. 531) 

0. 01397 
(1. Oil) 

0. 01257 
(1. 085) 

0. 0085 
(0. 180) 

0. 0073 
(0. 209) 

0. 0055 
(0. 398) 

0. 00495 
(0. 427) 

0. 01334 (0. 00525) 

500 x 500 

0. 001 

0. 00277 
(1. 090) 

0. 00254 
(1. 643) 



0. 00109 
(0. 429) 

0. 001 
(0. 647) 



0. 00254 (0. 001) 

150 x 150 

0. 0026 

0. 00838 
(0. 612) 

0. 00749 
(1. 278) 

0. 00772 
(1. 974) 

0. 00610 
(2. 639) 

0. 0033 
(0. 241) 

0. 00295 
(0. 503) 

0. 00304 
(0. 777) 

0. 0024 
(1. 039) 

0. 0066 (0. 0026) 

60 x 60 

0. 0075 

0.02743 
(0. 587) 

0. 02045 
(1. 252) 

0. 02997 
(1. 948) 

0. 02045 
(2. 614) 

0. 01080 
(0. 231) 

0. 00805 
(0. 493) 

0.0118 
(0. 767) 

0. 00806 
(1. 029) 

0. 01905 (0. 0075) 

40 x 40 

0. 01 

0. 0396 
(0. 455) 

0. 0368 
(0. 554) 

0. 0290 
(1. 034) 

0. 0247 
(1. 133) 

0.0156 
(0. 179) 

0. 0145 
(0. 218) 

0. 0114 
(0. 407) 

0. 00974 
(0. 446) 

0. 0254 (0. 01) 

a Reference 

20. 










b Values in 

parenthesis are channel spacing 

m cm. 








c Values in 

parenthesis are channel spacing 

tn inche s. 














of cogent reasons for the data scatter. First, it was extremely difficult to 
obtain accurate data using as the flow medium a near-saturated liquid, 
which must undergo a phase change prior to use in the measurement system, 
i. e. , the LH 2 in the apparatus must be converted to warm GH 2 outside the 
apparatus before pressurizing the manometers. In many cases, pulsing of 
the pressure was noticed, especially at the lower flowrates. In general, 
most frictional pressure drop experimentation is done with ambient- 
temperature low-pressure gases, to obviate this problem. Note that the low- 
pressure GN 2 data, at about R = 5, 000, show very little scatter relative to 
the LH 2 data. Another severe problem was the use of a high-pressure flow 
medium. Using LHg at 34. 5 N/cm^ (50 psia) greatly magnifies the effects of 
very minor leakage when determining pressure drops of perhaps 0. 0007 N/ 
cm^ (0. 001 psi). Extensive effort was expended during the test program to 
find and eliminate leakage from all sources, including the instrumentation. 

Substitution of equation (12) into equation (10) and examination of the 
Moody graph reveals another potential source of error. The friction factor 
is a very weak function of the relative roughness, where e may be rather 
accurately specified (for screens), but varies essentially as the cube of the 
channel spacing, s, which may be neither well-defined nor constant along the 
channel, due to normal manufacturing tolerances. Further, in the turbulent 
regime, the friction factor varies as the square of the flowrate, which showed 
variations of a few percent during our tests. All of these considerations 
may contribute to the observed data scatter. 

Another interesting feature observed from Figure 21 (and Appendix B) 
is that the data for the small channel spacing (s < 0. 64 cm (0. 25 in. ) - see 
Table IX) appears high compared to the computed e/D. It is thought that 
this may be due to the large w/s ratio for these cases. In laminar flow, 
solution of the equations of motion leads to equation (13) where the value of 
the constant, C, depends on the boundary conditions and the geometry of the 
duct. Eckert and Irvine (ref. 22) show curves giving the value of C for rec- 
tangles, triangles, ellipses, and annuli. While for circular ducts, C = 64, 
for very elongated rectangles or for thin annuli, C approaches 96, assuming 
f and R based on hydraulic diameter. For rectangles increasing from 
w/s = 5 to W/s = 20 (the range of our conditions), the value of C increases 
from 7 6 to 90. Schlicting (ref. 23), in describing the work of L. Schiller 
and J. Nikuradse on turbulent flow in noncircular ducts, indicates that use 
of hvdraulic diameter for all of these ducts leads to accurate correlation in 
the turbulent regime. The Moody graph is for circular pipes and implies 
applicability to noncircular ducts when correlated with hydraulic diameter 
only. However, all normal noncircular ducts have values of C close to that 
for circular pipes, i. e. , a w/s = 3 rectangle has C = 68, equilateral and 
isosceles, triangles have C =* 53, etc. No data are known at this time for 
very flat rectangles or thin annuli in turbulent flow, however, it is con- 
jectured that the variation in f due to thin annuli or very flat rectangles which 
occurs in laminar flow may also persist into the turbulent regime and may 
account for the small upward deviation of our data from the Moody graph. 


Analytical Models for the Experimental Data 

The correlation of f and R, with e based on screen characteristics, 
described previously was based on flow data from rectangular channels. To 
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use the data correlations in the analysis of the wall screen liner, which is a 
very thin annulus with no sidewalls (only a screen side and a tank wall side), 
either the data must be corrected for the effects of the sidewalls, or the 
effects of the sidewalls must be shown to be insignificant. An analysis to 
verify that the sidewall contribution to the total channel pressure drop can 
be ignored is shown in Appendix C. The results of this analysis indicate 
that at worst, the sidewall effect is much less than the data scatter, and at 
best, is insignificant. 

In addition, while our data correlate well with the Moody graph in the 
transition/turbulent regimes, the range of R for the system analysis task is 
from about 3 to over 100, 000. Thus, our correlation(s) must cover the 
complete range of flows in the laminar, critical, transition, and turbulent 
regimes . 

In the laminar regime, the correlation should clearly be, or approxi- 
mate, equation (13) with C = 96 for thin annuli (the case for our tankage sys- 
tems). In the critical regime, no correlation is available, and in the 
transition/turbulent regimes, our data are correlated by the Colebrook 
function of equation (15). 

It can be seen that at large R, the f values approach the Von Karman 
rough pipe formula: 


J. 


2 log 


3. 7 

e/D h 


(16) 


Determination of the friction factor from equation (15) is awkward, since 
the equation is only implicit in f. Further, in our tankage system analysis, 
for steady flow around a spherical annulus of constant width, the velocity in 
the annulus will be varying with angular position in the annulus. Therefore, 
both f and R will be varying with velocity, and equation (15) should be inte- 
grated around the spherical annular flow field. While this integration could 
be performed numerically, it would be very time consuming for the many 
cases of our extensive analysis matrix. In addition, since channel flow loss 
is only part of the total pressure loss, which includes head, velocity, and 
screen flow-through losses, such analytical tedium is unwarranted. Also, 
solution of equation (15) would still not solve the problem of the determination 
of f in the critical regime. 


Therefore, a similar correlation to the screen flow-through loss corre- 
lation was assumed for the channel flow loss. This correlation was asymp- 
totic to the laminar function, equation (13) at low R and to the rough pipe 
formula equation (16), at high R, and is: 



(17) 


This expression is shown in Figure 23 compared to the laminar function, 
equation (13), and to the Colebrook function, equation (15). At low values of 
e/D^, typical of most of our higher R cases, equation (17) is an excellent 
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approximation of equation (15), gives conservative f values in the critical 
regime, and closely approximates equation (13) at R < 103. At very coarse 
roughness values (e/Dh = 0.05), the f values are more conservative, but in 
our analysis matrix, the coarse screens are usually appropriate to systems 
with very low R, where the deviation from the laminar expression is small. 


Another consideration merits discussion. In ordinary pipe flow, where 
C = 64, the critical jump from laminar to turbulent flow occurs usually at an 
R of 2, 000 to 4, 000, and since f is small (0. 032-0. 016) the jump in f is 
upward to the turbulent value. However, the critical regime can start at R 
as low as 1, 000, depending on geometry, initial flow disturbances, entrance 
length, flow bends, etc. (References 20 and 22). It is conjectured here that 
the flow through the screen at the beginning of the channel flow annulus, 
though small, could act as a turbulence generator, thus triggering transition 
to turbulent flow at fairly low R ( < 10^), At this low R, f would have a value 
of 0. 096 or above, which is above the f values for turbulent flow. It is not 
reasonable that at transition the f values would drop abruptly to the values 
predicted by the Colebrook function, equation (15), but rather that the f 
values would make a smooth transition to the turbulent regime. 

The correlation of equation (17) provided this kind of smooth transition 
from laminar to turbulent flow, and could also be expressed with pressure 
drop as an explicit function of velocity in the form: 


AP 


96p 

2g D, 2 
6 c h 



(18) 


which is a convenient form for the Task 3 Analysis. 
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ANALYTICAL EVALUATION OF TANKAGE SYSTEMS 


Tankage System Operation and Weight Parameterization 

Prior to the analytical study of the tankage system, the operational 
aspects of the system were assumed as follows: during LIT? inflow to the 
empty tank, an auxiliary fill vent valve operates while filling the annulus 
and standpipe in low-gravity ( 10 - 5 g's with the g-vector toward the inflow 
line), also cooling all tankage and internal structure (Figure 1). The fill 
vent will be closed after the annulus and standpipe are full, and inflow will 
continue, nonvented, with pressure rise to 34. 5 N/crn^ (50 psia). During 
storage for periods of 30 and 300 days, the TVS is assumed to operate 
continuously, at a g- level of 10“ 5 g's with the g-vector toward the inflow 
line, and it is assumed that the tank is an adiabatic system at a temperature 
of 25. 2 0 K (45. 4 0 R) with suitable control to maintain a constant tank pres- 
sure of 34. 5 N/cm^ (50 psia) with hydrogen liquid and vapor as the only 
contained fluids. During outflow at constant pressure, the g-level vector of 
10“ 5 g' s i s away from the outflow line. 

All six tanks shown in Table I were studied at the following specified 
flowrates: 

• Inflow Rate - 1% (of tank volume) /min 

• Outflow Rate - 1 and 0. 01% (of tank volume) /min 

• TVS Flowrate - 1 and 0. 1% (of tank volume) /min 

The entire tankage system shown in Figure 1 was examined to determine 
the weight sensitivity of various system components for a given tank. The 
components and their functional weight relationship are shown in Table X. 
There are a number of discrete categories of related weight components. 

For example, items 5, 7, and 13 are only minor functions of things other 
than tank size, and thus can be considered fixed for a given tankage system 
and will not affect the weight optimization. The annulus residual and puddle 
residual are functions of the standoff distance (annulus gap) and screen type, 
as is the weight of the wall screen liner; these are parameterized and 
optimized in the next section. Items 10 and 11, the vent loss (boiloff) from 
external heat leak and the tank insulation weight, are interrelated and inde- 
pendent from the optimization analysis except as the total vent (boiloff) rate, 
together with the specified TVS flowrate (1 %/minute or 0. 1 %/minute), affects 
the TVS heat exchanger weight. The weight analysis for these items will be 
described in a later section. 

The items remaining (2, 6, 8, and 9) are all functions of pump size and 
will influence the overall system weight optimization. These items will be 
optimized in the section on determination of pump power requirements. 

The final section in the analysis will use the results of the above optimi- 
zation analyses to determine the optimum system configuration in terms of 
minimum weight. 
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TABLE X. - SYSTEM WEIGHT FUNCTIONAL RELATIONSHIPS FOR A GIVEN TANK 


Weight component 

Weight a major function of 

Weight a minor function of 

1. Annulus residual 

Standoff distance 


2. Standpipe residual 

Pump size 


3. Puddle residual 

Standoff distance, screen 


4. Wall screen liner 

Screen 


5. Wall screen liner supports 


Screen 

6. Standpipe 

Pump size, material 


7. Standpipe supports 


Standpipe size 

8. TVS pump/motor 

Pump size 


9. Boiloff from TVS pump 

Pump size, mission time 


10. Boiloff from external heat leak 

Mission time, insulation type 


11. Insulation 

Mission time, insulation type 


12. TVS heat exchanger 

TVS pump flowrate 

Standpipe size, boiloff rate 

13. Baffles 


Tank diameter 




Screen Size and Wall Spacing Analysis 

To evalute the performance of the wall screen liner system for LH 2 
acquisition during outflow, the annulus flow correlation and model developed 
from the channel flow loss data in the previous section was adapted to a 
spherical annulus, such as shown in Figure 24. Equation (18) described 
pressure loss, AP, as a function of velocity V, (see symbols): 


AP = 


96 p 

2g D? 
s c h 


LV + 


LV 


8 (log 


3. 7 
e/D. 


g D. 
B c h 


Equation (18) was integrated around the spherical annulus, shown in 
Figure 24 where 

Dt 

L = ~2 d^ (19) 

in polar coordinates. 

Equation (18) became: 


AP 


A 



D r 


V 


d 0 + 





V 



( 20 ) 



F igure24. Polar Coordinate System for Spherical Screen Annulus Model During Outflow 
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where 


A = 


96 |j. 

? * 

2 g D i 

°c h 


B = 


8 < lc, g fyD> S c D h 


For steady volumetric flow. 


Q - VsirD^cos j> - constant 


( 21 ) 


Substituting in equation (20) gave: 


AP 


AQ 


2 TT S 


i- / 

G J 


01 


2 BQ 

cos 0d0 + 2 — 2 


2 tt s D„ (f . i 


/ 2 2 
J COS 0d0 


(22) 


and integrating gave: 


AP = 


96 [j. Q 


2 g c D h 2TTS 


(^sin + sin + 


PQ 


2. 2 s 2 D T 8(log^-) g c D h 

h 


0^ s i- n 2^2 s in 20^ 
— - + ; 


2 + -f + 


(23) 


This equation was, in fact, erroneous, because it was based on rectilinear 
flow in a cylindrical annulus which was artificially integrated in spherical 
coordinates in a way which did not properly account for the curvature of the 
streamlines in the spherical annulus. To obtain the correct expression, the 
equations of motion must be solved in spherical coordinates with the appro- 
priate boundary conditions. This solution can only be obtained in closed 
form for very slow laminar flow, i. e. , such that the inertial terms in the 
momentum equation can be neglected. This analysis is shown in Appendix D, 
and, dimensionally correcting jj. in Appendix D to p/g c and using the nomen- 
clature of Figure 24, gives the expression: 

f A. r <f> 7 + tt/2 / u/2 -0,-1 

AP = ^ — -J In Ttan 2 j tan 2 (24) 

g c it s L / -* 


The first term, or laminar portion, of equation (23), with = 2 s for an 
annulus, is 


AP = ^ [sin ^2 + S ^ n ^ll 


(25) 


g TT S 
®C 
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The difference in the two equations is confined to the final trigonometric 
terms, which, for comparison, are plotted with = <f>\ in Figure 25. The 
functions are essentially equal at angles less than tt/6 radians, but diverge 
widely at angles near tt/2 radians. At angles of 5 ir/12 radians, which are 
typical values for our spacing analysis discussed below, the correct value 
is twice the previously modeled value. While the correct function could be 
substituted in the laminar component of equation (23), no equivalent expres- 
sion exists for the turbulent component of equation (2 3). 

As an approximation, it was assumed that the first part of equation (23) 
accounts for viscous flow (that which ignores inertia terms in the momentum 
equation) while the second part of equation (23) accounts for the inertial 
effects. At the same time the ratio 



corrects for the geometric conditions of curvature of the streamlines in 
laminar flow. Streamlines would also be present in turbulent flow, even 
though velocity- dominated turbulence would be superimposed on the stream- 
lines. Curvature of the turbulent streamlines because of geometry should 
require the same kind of pressure correction. Therefore, it was assumed 
that the entire pressure loss described in equation (23) be multiplied by the 
geometry ratio of equation (26). This would give a reasonable approxima- 
tion of the pressure loss in the spherical annulus. In the straight annulus 
of tanks with an L/D ratio above 1, the original correlation for a cylindrical 
annulus (without the above correction) was used. It should also be noted 
that the spherical annulus flow loss is only part of the total flow loss which 
includes static and dynamic head loss and loss through the screen. Thus, 
errors in the annulus flow loss would be less significant relative to the total 
loss. Further, the screen spacing analysis described below indicated that 
screen flow-through loss, rather than annulus loss, was strongly dominant 
for most configurations studied. 

Expressions for static head loss, dynamic head loss, and residual were 
also derived. Referring to Figure 24, the configuration shown is the worst 
case, since the static head is maximum, the channel length is maximum, 
and the screen flow-through area is minimum. This situation would occur 
at the end of outflow, with the tank nearly empty, and screen breakdown 
about to occur at the outflow baffle. The TVS flow case would be exactly 
reversed, the g-vector reversed, and the flow going from the outflow baffle 
to the other baffle. Because screen flow-through loss does not occur, the 
TVS case is not as severe as the outflow case. 

The g-level during outflow was specified as 10“^ g's. The static head 
loss, Hg, was thus the g-level times the length from outflow baffle to liquid 
surface, or: 


H c 


_ 5 

10 

- 1 

) D T + 

(sin + sin <^> 2 ) 

(D T /2 - s) 







(27) 
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Since breakdown would occur at the screen next to the outlet baffle, the 
dynamic head loss, H^, at the outlet baffle was determined by the velocity 
in the annulus at the outlet baffle location, or: 


H. 


V 2 /2g, 


where 


V = 


Q 

7tD 2 S 


(28) 


The pressure loss through the screen, H s , was determined from the flow 
loss correlation described in equation (7) and Table VIII. 

H s = AVj + BVj 2 (29) 

where A and B are the experimentally determined coefficients from Table VIII, 
and Vj is the velocity through the screen, which depends on the flowrate and 
the spherical segment screen area. A, where Vj = Q/A and: 

A = 2 tt (sin T - sin <f> (D T /2 - s ) 2 (30) 

The head loss in the annulus was added to the screen flow-through loss, 
the dynamic head loss, and the static head loss to give the total loss. The 
screen bubble point was divided by the total head loss to give the safety 
factor. 

The residual in the tank was made up of the annulus residual, the stand- 
pipe residual, and the puddle residual, i. e. , the residual in the tank bottom 
between the standpipe and annulus. The annulus residual, Xj, for a thin 
annulus, was fixed for a given annulus thickness, s, and is: 

x i = <5^> 5 ” d t 2 < 31 > 


The standpipe residual, X 2 , was fixed for a given standpipe diameter, 
Dj and is: 


X 


2 



D T - 2 s 



(32) 


The puddle residual varied with the angle and was the volume in the 
spherical segment between the standpipe and annulus. A flat interface was 
assumed for convenience, and gives conservative results, since with a curved 
interface with wetting LH 2 , the residual would probably be less, depending on 
the interface curvature. (In practice, for our analysis cases, the puddle 
residual was very small, as will be discussed below, so that this conserva- 
tism did not materially affect the study conclusions.) The puddle residual, 

X 3 , is: 
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Y 


(33) 


X 3 = T Y 


(3 + 


Y 2 ) 


ttD 


where 

Y = (D t /2 - s) (1 - sin 0 ) 

and 

Z = (D t /2 - s) (cos 0^) 


The computations were organized as follows: first, the fluid properties 
of density and viscosity were input, then the tank parameters-volume, L/D^, 
standpipe, and baffle diameters were input. Next, the screen characteristics 
of bubble point, flow-through coefficients (A and B in equation (29)), and 
roughness dimension, e, were input. Next, the annulus gap (screen standoff 
distance) and flowrate were input. The program computed and printed the 
baffle angle, T, and asked for an input of which had to be less than T to 
give a positive flow-through area (Figure 24) plus an input flag to instruct 
the program as to whether the next case would be with a new angle 0^, a new 
annulus gap, a new screen, or a new tank, baffle, and standpipe configuration. 
The program then computed and printed tank diameter, submerged screen 
length, standpipe residual, annulus residual, puddle residual, total residual, 
percent residual, static head loss, dynamic head loss, screen flow-through 
head loss, annulus head loss, total head loss, and safety factor. The program 
was coded for use on the MDAC Time-sharing Computer System and was 
straightforward in operation: for a given tankage, screen, annulus gap, and 
flowrate combination, a variation in gave a variation in puddle residual 
and safety factor which could then be crossplotted to determine system 
sensitivities. 


The six tanks selected for the analysis were listed in Table I, and will 
be identified subsequently by volume over L/D ratio, i. e. , 5, 000/4. Each 
tankage system was analyzed to determine the choice of standpipe and baffle 
diameters. Since the standpipe residual is a fixed value which does not 
affect the screen wall spacing analysis sensitivities, the standpipe sizing 
was ignored, since it will be optimized along with the pump in the next sec- 
tion. The baffles at each end of the standpipe (Figure 1) were arbitrarily 
sized at about 1/4 the tank diameter as shown in Table XI. Since the outflow 
baffle is also part of the tank manhole for TVS system access. Table XI 
indicates that the arbitrary sizing results in reasonable baffle /manhole 
sizes and in very small minimum puddle residuals due to baffle size. 

Each tank was analyzed with all 10 screens and with annulus gaps from 
1% to 5% of tank volume, as shown in Appendix E. The screen performance 
parameters are summarized in Table XII. Typical results of this analysis, 
for the 5, 000/4 tank, are shown in Figures 26 and 27. Figure 26 shows the 
performance of the 10 screens in terms of safety factor and puddle residual 
for the 1% annulus gap (the specified minimum). Since the avowed purpose 
of the ultimate tradeoff study was to minimize system weight, it was logical 
to select screens with adequate performance at the minimum annulus gap. 
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TABLE XI. - TANKAGE SYSTEM AND BAFFLE PARAMETERS 


Tank ID 

Tank diameter, 
m (ft) 

Upper and lower 
baffle diameter, 
m (ft) 

Minimum puddle 
residual, 

% 

5,000/4 

3. 66 (12. 02) 

0. 91 (3. 0) 

0. 0104 

500/4 

1. 70 (5. 58) 

0. 46 (1. 5) 

0. 014 

500/2 

2.21 (7.25) 

0. 55 (1. 8) 

0. 022 

500/1 

3. 00 (9. 85) 

0.73 (2.4) 

0. 052 

50/2 

1. 03 (3. 37) 

0. 27 (0. 9) 

0. 03 

50/1 

1. 39 (4. 57) 

0. 37 (1.2) 

0. 07 


For system/screen comparison, a safety factor of 2 was chosen as 
representative of adequate performance. Figure 26 shows that at a 1% annu- 
lus gap, only the 325 x 2, 300 and 200 x 1, 400 screens gave a safety factor 
of 2. Since the 325 x 2, 300 screen was substantially lighter than the 
200 x 1, 400 screen, and had a higher bubble point (Table XII) it was the 
logical choice. Note that the puddle residuals obtained were very small and 
in fact insignificant compared to the 1.416 m^ (50 ft^) residual in the 1% 
annulus . 

From Table XII, there are only two screens lighter in weight than the 
325 x 2, 300 screen: the 500 x 500 and 150 x 150 screens. The 500 x 500 
screen is extremely flimsy, while the 150 x 150 screen appears to have 
reasonable structural rigidity, perhaps comparable to the 325 x 2, 300 
screen. Figure 27 shows that at 2% annulus gap, the 150 x 150 screen had 
comparable performance to the 325 x 2, 300 screen in the 1% annulus. Thus, 
for the 5, 000/4 tank with 168. 5 (1, 814 ft^) of s creen surface area, the 

150 x 150 screen would save 32 kg (71 lb) of screen weight, but at a cost of 
1% or 91 kg (200 lb) in annulus LH 2 residual weight — a net increase of 
58. 5 kg (129 lb) in total weight — assuming comparable structural qualities. 
The 500 x 500 screen, in a 1.20% annulus, would have comparable perform- 
ance to the 325 x 2, 300 screen in a 1% annulus, and the 18 kg (40 lb) 
increase in annulus LH 2 residual weight would be more than offset by the 
62 kg (136 lb) saving in screen weight — for a net decrease of 44 kg (96 lb) 
in system weight. However, the extreme flimsiness of the 500 x 500 screen 
would certainly require some sort of increased structural backup material 
for installation in a flight vehicle, and the weight of this increased backup 
material would certainly exceed the 44 kg (96 lb) of weight gained. Thus, 
for the 5, 000/4 tank, the 325 x 2, 300 screen would be the logical choice on 
the basis of weight and performance 

Appendix E indicates similar results for the 500/4 tank, where 5 screens 
have adequate performance; for the 500/2 tank, where 8 screens have ade- 
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TABLE XII. - SCREEN PERFORMANCE PARAMETERS 


Screen 

Bubble point, 
m (ft) LH 2 

Flow-through parameters 

Roughness, 

Weight, 
kg/m 2 
(lb/100 ft 2 ) 


A 

B (ft) 

€ 

cm (in. ) 

325 x 2, 300 

0. 4815 

1. 14 

0. 6919 

0. 00127 

0. 532 


(1. 580) 


(2.27) 

(0. 0005) 

(10.9) 

200 x 1,400 

0. 3377 

0. 885 

0. 6126 

0. 00203 

0. 908 


(1. 108) 


(2. 01) 

(0. 0008) 

(18. 6) 

720 x 140 

0. 1767 

0. 162 

0. 2627 

0. 00546 

0. 693 


(0. 580) 


(0. 862) 

(0. 00215) 

(14.2) 

165 x 800 

0. 1228 

0. 108 

0. 0805 

0. 00254 

0. 796 


(0. 403) 


(0. 2 64) 

(0. 001) 

(16.3) 

50 x 250 

0. 0682 

0. 045 

0. 0631 

0. 00572 

1. 133 


(0.224) 


(0.207) 

(0.00225) 

(23. 2) 

24 x 110 . 

0. 0336 

0. 0165 

0. 1554 

0. 01334 

3. 1 


(0. 1105) 


(0. 51) 

(0. 00525) 

(63. 5) 

500 x 500 

0. 1646 

0. 0554 

0. 04267 

0. 00254 

0. 166 


(0. 540) 


(0. 140) 

(0. 001) 

(3.4) 

150 x 150 

0. 0460 

0. 0101 

0. 01347 

0. 0066 

0. 342 


(0. 151) 


(0. 0442) 

(0. 0026) 

(7.0) 

60 x 60 

0.0230 

0. 00585 

0. 02761 

0. 01905 

1. 167 


(0. 0754) 


(0. 0906) 

(0. 0075) 

(23. 9) 

40 x 40 

0.0170 

0. 00287 

0. 01497 

0. 0254 

1. 362 


(0. 0559) 


(0. 0491) 

(0. 01) 

(27.9) 
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quate performance; while all 10 screens would have adequate performance 
for the remaining tanks. The weight savings using the 150 x 150 screen 
instead of the 325 x 2, 300 screen ranged from 5. 9 kg to 1. 1 kg (12. 9 lb to 
2. 5 lb) (which were not significant weight savings from a system standpoint) 
while sacrificing an order of magnitude in bubble point. 

The analysis described above was for the maximum outflow rate of 1% of 
tank volume /minute. At 0.01% outflow rate (or at 0. 1% TVS flowrate) there 
was absolutely no sensitivity to screen type or annulus gap; all 10 screens 
had adequate performance at the minimum 1% annulus gap. It is probably 
not possible to significantly reduce the annulus gap below 1%, since the 
1% annulus varies from 0. 84 cm (0. 33-in. ) in the 141. 6-m^ (5, 000-ft^) tank 
to about 0.23 cm (0.09 in.) in the 1.416-m^ (50-ft^) tank; these are small 
values from a system fabrication standpoint. It is clear that at this g-level, 
flowrate, and annulus gap, the tankage systems flow characteristics were 
such that a meaningful tradeoff analysis in terms of residual, annulus gap, 
and performance was not possible. 

Therefore the analysis was extended to examine the effects of high 
outflow rate (3% tank volume /minute) and high g-levels (10"^ to 10-4 g's) on 
system performance. Only the high L/D tanks (5, 000/4, 500/4, and 50/2) 
were analyzed because these were generally the more severe cases. 

Figure 28 shows the effect of increasing the flowrate to 3%/minute on the 
performance of the 325 x 2, 300 screen in the 5, 000/4 tank. The annulus 
gap had to be increased to 2% to achieve adequate performance, and even 
then, increased residual resulted. The results for the other tanks are also 
shown in Appendix E. 

_5 

A similar trend occurred when the g-level was increased from 10 g's 
to 10"^ g's, as shown also in Appendix E. At 10“^ g's, the 325 x 2, 300 
screen no longer had adequate performance in the 1% gap in the 5, 000/4 tank, 
and the annulus gap had to be increased. In Figure 29, the annulus gap is 
2% in the 5, 000/4 tank. The 325 x 2, 300 screen had adequate performance 
at all g-levels, but the 150 x 150 screen did not have adequate performance 
at even 10"4 g's. If the annulus gap was increased to 3%, however, the 
150 x 150 screen would have adequate performance at even 10 “ 3 g's. At 
10-2 gi s> the hydrostatic head in the 5, 000/4 tank exceeded the bubble point 
of the 150 x 150 screen by a factor of 3, so that it could not practically be 
used. 

The conclusion reached from the study of increased flowrate and 
g-levels is that it is even more important that the lightest high performance 
(high bubble-point) screen (32 5 x 2, 300) be used. While the lighter 
150 x 150 screen may have adequate performance for minor deviations 
from the original design conditions, it does not have the bubble-point 
reserve capacity to handle extremes of flowrate or g-level with adequate 
performance. 


Determination of Pump Power Requirements 

The standpipe together with the annulus, remains full of LH2 following 
outflow. This LH 2 represents a residual weight penalty which it would be 
desirable to minimize. However, it was noted that arbitrary reduction in 
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Figure 28. High Outflow Rate Performance of 325 x 2,300 Screen in the 5,000/4 Tank 
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Figure 29. High G-Level Performance of Both Screens in the 5,000/4 Tank 
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standpipe size, which would decrease the weight of the standpipe and its 
residual, would also increase the pressure drop down the standpipe, which 
in turn would increase the power requirements of the TVS pump. This 
increased power requirement would lead to an increase in TVS pump/motor 
weight and in increased "boiloff" from the TVS pump/motor inefficiency. 

Therefore, it was clear that there was an optimum standpipe diameter 
which would minimize the combined weight of standpipe, standpipe residual, 
and pump boiloff due to pressure loss in the standpipe. The pump boiloff 
due to pressure loss around the annulus was not directly dependent on the 
standpipe diameter, did not enter this optimization, and will be accounted 
for later in the analysis. Similarly the pump/motor weight was a very 
small value, so that it too was ignored in the optimization, and will be 
accounted for later. 

The weight of the standpipe residual, in terms of the standpipe 
diameter, D s , and length, L, is (see symbols): 

irD 2 L 

W l = 1 P (34) 


The weight of the standpipe depends on the thickness of the standpipe 
and the material. Since there is essentially no pressure load on the stand- 
pipe, the thickness criterion used was that specified by NASA-MSFC as 
minimum handling gage for ducting in the Space Shuttle. The thickness in 
meters (inches) is: 


t AirT , T = 0. 000558 + 0. 024 D 
MIN s 

(tnyr-nvr = 0* 022 + 0. 024 D ) for Stainless Steel 
MIN s 

= 0* 00076 + 0. 036 D 
MIN s 

(t. ,,, T = 0.030 + 0.036 D ) for Aluminum 
MIN s 


(35) 


Multiplying the thickness by the density of Steel and Aluminum gives, 
for standpipe weight: 


W., = TT D L (A + BD ) 
2 s s 


(36) 


where 



A = 

4. 41, 

B 

A = 

2. 1, 

B 


190 for Stainless Steel 
99. 5 for Aluminum 


The weight of the "boiloff" due to the pump power is equal to the power 
dissipated in the LH;? times the mission time and divided by the LH 2 heat of 
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vaporization. It was straightforward to show that essentially all of the input 
power to the pump/motor is dissipated to the LH 2 , causing "boiloff" (see 
Appendix F). Therefore, the "boiloff" weight is 

w = _£L£JLt (371 

W 3 nh. J ( n 

1 fg 


where t is the mission time, h.£_ is the heat of vaporization, J is the energy 
conversion, 0. 102 kg-m/joule, and rj is the overall efficiency. The fluid 
power is QpH, where Q is the volumetric flowrate and H is the pressure 
drop (in m of LH 2 ) down the standpipe. This head loss is: 


H = f 


L V 

D 2g 
s s c 


In terms of the volume flowrate, Q = VA, or 

Q • 4 




TT D 


(38) 


(39) 


The friction factor, f, is a function of Reynolds number, R. For our 
flow conditions, the flow is turbulent and the standpipe hydraulically smooth 
so that the correlation of Blasius (ref. 24) is suitable, or: 


or since 


0.316 
“ R 0 - 25 


R 


4 p Q 

TT (JL D * 

n s 


0. 316 D 


0. 25 


f = 


(4 p Q/irp.) 


0. 25 


(40) 


(41) 


As shown in Figure 30, the Blasius correlation is accurate to within 
5% for R from 3, 000 to 300, 000. 


Combining equations (38), (39), and (41) gives: 


H 


0. 316 L (6 ' 4 hY 


(4 p Q/^hl) 0, 25 2g ™ 4,75 


D 
c s 


(42) 


W 


The total weight, from equations (34), (36), (37), and (42) is: 
2 


•rrD “ L 

W = 7 p + ttD L (A + BD ) 

r s s 


4 

Q pt 


0. 316 L (Q • 4/tt ) 
fg J \(4 p Q/rrp.)°‘ 25 2g 


1 


D 


4. 75 


(43) 
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Differentiating with respect to D and equating to zero gives: 

~ s 

2tt L {- + b) D 6,75 + ttLA D 5,75 
' 4 / s s 


-4. 75 


Qpt (0. 316) L (Q • 4/tt) 2 
n h fg J (4pQ/TTfj.)°* 25 2 g c 


0 


( 44 ) 


All of the parameters of equation (44) are known except D g and q. 
Equation (44) was solved for D s using the Newton- Raphson iteration tech- 
nique for various values of q and the results are shown for the six tanks and 
two values of TVS flowrate (1%/minute and 0. 1%/minute) in Appendix E, and 
typically as the nearly horizontal straight lines in Figure 31. 


To determine the value of D g which minimized system weight, the 
correct value of q was determined. This value of q was determined from 
the total pump power, not just that required to overcome the pressure loss 
in the standpipe. To circulate the TVS flow, the pump must provide dynamic 
head, hydrostatic head, and enough head to overcome frictional losses in 
the screen annulus and the standpipe. Actually, the hydrostatic head is 
recovered and converted to overcoming friction loss in the standpipe, but 
for conservatism this effect was ignored since the hydrostatic head in these 
cases was very small. Therefore, the total fluid power (Pf) required is 


P. = Q p (H + H + H , + H ) 
f K s g d a 


(45) 


To determine the overall efficiency, knowing the fluid power, the pump 
was thoroughly evaluated. The basis of this evaluation was the definitive 
study of cryogenic mixers by Poth, et al (ref. 2 5). 

One of the basic characteristic parameters associated with pump design 
and applications is the specific speed. The specific speed, N , is defined 
in this study as 

N = (n, r P m ) (G gpm) 1 ^ (46) 

S (H, feet) ' 


The specific speed as defined is dimensionless but an inconsistent set of 
units is used. 


A constant to account for unit conversion yields a consistent definition 
of the specific speed, N g ; the specific speed thus defined is 

N 1 = N /2815 (47) 

s s 

where 


_ (£7, rad/sec) (Q, m /sec) 

s / 2 i ZTTfe 

(g c H, m /sec ) 


1/2 
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N g was used when analytical derivations were required. The results 
are presented in terms of N s by the use of the appropriate conversion 
constant, 2,815. 

The specific speed characterizes the type of pump required, as shown 
in Figure 32. For our applications, requiring relatively large flowrate and 
low pressure rise, the axial flow pump with a high specific speed (> 8 , 000 ) 
was the proper choice. 

The head coefficient, ijj, is a dimensionless pressure (head) rise of 
the pump. The head coefficient, vp , is defined as 

+ = —TT <«> 

n r B 

where r^ is the blade radius. 

The pump pressure rise as a function of flowrate is a convenient way 
to present pump output performance. In nondimensional form, the pressure 
rise across the pump as a function of flowrate can be represented in terms 
of the head coefficient i|j and the flow coefficient, <f>. The flow coefficient, 

0 , can be defined in terms of the specific speed, the head coefficient and 
the geometry coefficient, £, as 

0 = (N g /2815) 2 ^ 3 / 2 /tt | 2 (49) 


where 

£ 2 = A 2^ TTr B 


and A 2 is the axial flow cross-sectional area of the pump. A vane- axial 
pump usually has a single rotor stage with a de swirl- stationary stage. The 
flow into the pump is usually assumed to have no prerotation. As a result, 
an approximate relation exists between the head coefficient and the flow 
coefficient. 


ip = 1.0-0 cot P 2 (50) 

where (3 2 the r °tor stage exit flow angle. Eliminating 0 from the above 
equation, a relation between N g and 0 is obtained so that 

N = 2815 (ir/cot 
s z 


1/2 


(r h /r B>‘ 




1/2 


(51) 




where r, is the hub radius, 
h 
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Typical values for the rotor stage exit angle, Pg = 0* 366 radian (21°), 
and for r^/rg = 0. 7, so that: 


4 ; = 


n 2/3 


4. 88 x 10 


(l-*!*) 


(52) 


Also, from equation (48): 



(53) 


With these equations, the characteristics of the pump were determined. 
There were two ways chosen to parameterize the performance of the pump 
for this study. One was to fix the specific speed at the maximum value 
possible for efficient vane-axial fans, 16, 000 (ref. 26), and allow the RPM 
to vary. With an AC motor, the RPM could be varied by altering AC fre- 
quency or poles in the motor or by gearing down the motor speed. The 
advantage of this method was that the fluid efficiency of the pump could be 
assumed because the specific speed and Reynolds number of the pump are at 
normal values. The efficiencies of various pumps are shown in Figure 32. 
For axial pumps, the efficiency shown is 92% at a Reynolds number of 106 , 
however, it is highly probable that the efficiencies shown are for very large 
machines. As pumps get smaller, fluid efficiency invariably suffers. 
Therefore it was assumed, for the very small vane- axial pumps required 
in this study that the maximum fluid efficiency was 70% at a specific speed 
of 16, 000 (ref. 26). The disadvantage of fixing specific speed would be 
that for low head requirements, the rotational speed was of necessity low, 
leading to larger pump diameter and heavier pump weight, however, the 
pump diameters turned out to be very close to the standpipe diameters, 
which was desirable from a system integration standpoint. 

The other pump parameterization method was to fix the pump rotational 
speed at a value suitable for a normal electric motor, say 100 v radians/ 
sec (3, 000 rpm) and allow the specific speed to vary. Unfortunately, for 
low head rise machines the resulting specific speed was very high, approach- 
ing 105 to 10&, and the performance of such high specific speed machines is 
invariably severely degraded. In fact, the efficiency of axial pumps drops 
rapidly at specific speeds above 20, 000 (ref. 26). Therefore, because 
use of this pump parameterization method resulted in uncertain (but very 
low) efficiency definition, it was not used in the final optimization analysis. 

Estimates of the overall electric pump motor efficiency as a function 
of fluid power were required to determine the weight attributable to the 
pump/motor subsystem. Data were obtained from Stark (ref. 3) and 
Sterbentz (ref. 2) in previous studies conducted in this area. The results 
are shown in Figure 33 for both AC and brushless DC motors. 

Based on the data available, the DC brushless motors have a higher 
efficiency, especially at low fluid power. No actual data on the brushless 
DC motor operating in liquid hydrogen were available, whereas data for 
AC motor efficiency were obtained from actual LH^ pumps (ref. 26). As a 
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Figure 33. Comparison of Electric Motor Efficiencies 


result, the AC motor data were used in the study. A curve fit of the AC 
motor data at low power resulted in the following equation: 

0 524 

r\ = 0. 121 (P.) (54) 

where q e is the electric motor efficiency, and P| is the input power in watts 
to the motor. Assuming a fluid efficiency of 70%, the overall efficiency, 
r|£, in terms of the fluid power in watts, P£, is: 

0 544 

n t = 0. 199 (P f ) “ (55) 

The correlations for pump and motor weight were also taken from curve fits 
obtained from References 2 5 and 2 6. The pump weight correlation in kg 
includes a 50% factor for additional weight to integrate the pump with the 
outflow baffle, and is: 

W p = 245 D b 2 * 34 (56) 

the fan blade diameter in meters. 


motor weight (kg) is 


where is 
The electric 



where r ] £ is the fluid efficiency (set at 70%) and n is the motor rotational 
speed in rpm. 

From equations (42), (45), and (55), the overall efficiency as a function 
of the total fluid power, including the standpipe loss, was determined and 
plotted typically in Figure 3 1 and Appendix E. Where these lines crossed the 
standpipe optimization lines, the intersection was the value of standpipe 
diameter which gave the minimum system weight for that annulus gap, 
screen, and standpipe material. It was assumed that either stainless steel 
or aluminum could be used for the standpipe because of the presence of the 
slipjoint in the standpipe (Figure 1) which allowed differential thermal expan- 
sion. With the design value of standpipe diameter as an input, and with the 
other head losses known, the total system weight analysis was completed. 


System Weight Optimization Analysis 

The overall system weight was divided into four categories: annulus 

residual, standpipe residual, pump boiloff, and hardware, (which includes 
standpipe weight, screen weight, pump and motor weight. ) The puddle 
residual, as shown in the annulus gap analysis, was insignificant and was 
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ignored. The tankage insulation system and externally- caused "boiloff" 
were optimized for both the 30 and 300-day mission as shown in Appendix G. 
They had no effect on the optimization except as they affected the TVS heat 
exchanger weight, however, it is of interest to note that after 300 days, the 
5, 000/4 tank had lost 11% of the stored LH 2 to external boiloff, while the 
500/4, 500/2 and 500/1 tanks had lost about 22%, and the 50/2 and 50/1 tanks, 
about 43% to external boiloff. As shown in Appendix G, the largest heat 
exchanger, for the 5, 000/4 tank, consisted of about 6 m (20 ft) of 0. 65-cm 
(1/4-in. diameter) tubing and weighed less than 0. 7 kg (1. 5 lb); therefore, 
the effect of the TVS heat exchanger optimization and design on the overall 
system weight was ignored. The four categories above plus the total weight 
is shown for the six tanks, for 30 and 300-day missions, for various screens 
and standpipe materials versus annulus gap in Appendix E, as typically shown 
in Figures 34 and 35. The results for the 30-day mission are shown typi- 
cally in Figure 34 and indicated that minimum system weight was achieved 
with the minimum annulus gap because of the strong influence of annulus 
residual. 

However, for the 300-day mission, the pump boiloff became very 
important, as shown in Figure 35, so that an optimum annulus gap was found, 
at which minimum weight occurs. For the 5, 000/4 and 500/4 tanks, this 
optimum was at about 2. 0 to 2. 2% annulus gap for both the 325 x 2, 3 00 and 
the 150 x 150 screens. For the 500/2 and 50/2 tanks, the optimum gap was at 1 
about 1. 6%, and for the 500/1 and 50/1 tanks the optimum was at an annulus 
gap of about 1.4 to 1. 5%. For the 5, 000/4 tank at 0. 1% tank volume /minute 
TVS flow, the optimum again occurred at the minimum annulus gap, as was' 
the case for the other five tanks. 

The optimum aluminum standpipe resulted in a weight reduction of 8 to 
10% of the total compared to the optimum stainless steel standpipe, and use 
of the 150 x 150 screen resulted in a further weight reduction which was 
entirely due to use of a lighter screen material. 
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CONCLUSIONS 


The conclusions reached from the wall screen spacing analysis for out- 
flow are that the 325 x 2, 300 screen gave maximum performance and reserve 
safety factor compared to the other nine screens but that the lighter 150 x 150 
screen could be used with, at worst, a larger annulus gap. The 150 x 150 
screen showed reduced performance compared to the 325 x 2, 300 screen at 
higher flowrates and g-levels. However, the results of the weight optimi- 
zation based on the TVS pump operation for long storage times may alter the 
conclusions based on outflow. For example, with the 5,000/4 tank and 30-day 
storage time, the minimum would be achieved with the 325 x 2, 300 screen 
and 1% annulus gap, because the lighter 150 x 150 screen would not have 
adequate performance. If, however, the system were designed for a 300-day 
storage time, the minimum weight design would be with a 2% annulus gap, 
and the lighter 150 x 150 screen could be used without sacrificing outflow 
performance. 

In our study, both the analytical and experimental work were constrained 
to Type 304 stainless steel (test apparatus and model tank) so that the 
screens selected were also of stainless steel for material compatibility. 
However with aluminum tankage, a further major reduction in screen weight 
could be obtained by using aluminum screens. 

Further conclusions which may be drawn from these results are that for 
moderate missions (30 days or low TVS flowrate) where the pump boiloff is 
not really important, inaccuracies in the description of the pump/motor 
efficiency will not have a marked effect on the results; which, predictably, 
favor use of the minimum available annulus gap. On the other hand, for the 
long duration missions (300 days) the determination of pump efficiency is 
most important because pump boiloff becomes a dominant effect and drives 
the design of the TVS/WSL system. The curve fit used for the efficiency 
evaluation was extrapolated over several orders of magnitude from rather 
scattered data (Figure 33). For the 5, 000/4 tank, the input power level was 
several watts, and there can be some confidence that the efficiencies and 
power levels are reasonable, and the pump and motor designs achievable. 

For the other, smaller tanks, the power levels were so low that it is not 
clear that the pump/motors can be built or that the predicted efficiencies can 
be achieved. Further definition of practical pump/motor configurations, 
together with a more confident evaluation of efficiencies should be determined 
before design of the TVS/WSL system for a 300-day mission be undertaken. 

In addition, operational alternatives, such as operating a larger, more 
efficient pump/motor intermittently to reduce overall heat imput, should be 
evaluated. This would require definition of in -tank flow fields and thermo- 
dynamics in low gravity, and investigation of alternate ways of controlling 
tank heat input, such as vapor-cooled shields around the tank. 

The conclusions reached from the experimental phases are: 

A. The bubble points of the 10 representative screens in 34. 5 N/cm 4 
(50 psia) saturated LH 2 were accurately predicted from bubble- 


85 



point tests with isopropyl alcohol. The bubble-point data corre- 
lation was generally lower ( — 13%) than the correlation based on 
other LH 2 data. 

B. The screen flow-through pressure loss data taken with both 

ambient GNg and 34. 5 N/cm.2 (50 psia) saturated agreed well 

with other flow loss data, and the correlations obtained spanned the 
full practical range of Reynolds numbers. 

C. The channel flow-loss data taken with 34. 5 N/cm^ (50 psia) LH 2 
were unique, but the data were well correlated with the friction 
factor based on a roughness parameter determined from the screen 
shute wire diameters. 

A final conclusion is that for most tanks and mission times, the annulus 
residual and hardware weight (mostly consisting of the wall screen liner) 
account for the great majority of the system weight penalty. This is a direct 
consequence of using a full wall screen liner as the acquisition system. 

Now that this study has shown the overall concept to be fluid -dynamically 
feasible, especially for large tanks, it is recommended that problems of 
tank thermodynamics and thermal control be investigated, with the hoped-for 
result of obtaining realistic additional system optimization and weight 
reduction. 
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APPENDIX A 


SCREEN SURVEY TABULATION 
WIRE CLOTH TERMINOLOGY 


1. Bolting Cloth 

2. Calendered Cloth 

3. Count 

4. Dutch 

5. Market Grade 

6. Plain Dutch 

7. Reverse Plain Dutch 

8. Selvedge 

9. Shute 

10. Sintered Mesh 

11. Square Weave 

12. Twilled 

13. Warp 


A precision-woven square-weave cloth woven 
on special loops of custom drawn wire with a 
very smooth finish. 

Wire cloth which has been passed through a 
pair of heavy rollers to reduce the thickness, 
reduce pore size, or flatten the cloth. 

Number of openings per lineal inch. 

Warp and shute wires of different diameter. 

An economical straining square -weave cloth. 

Similar to square weave, except warp wires 
are heavier and the shute wires are driven 
closely together. 

Warp wires are placed close together and the 
heavier shute wires are woven tightly over one, 
under one. 

A finish edge on wire cloth to prevent ravelling. 

Wires running crosswise in the cloth as woven. 
They are passed back and forth through the 
warp wires by the shuttle of the loom. 

Woven wire mesh made more rigid by furnace 
bonding the wires of the weave at all contact 
points . 

Square mesh has warp and shute wires of 
equal diameter spaced equally in both direction 
directions. Wires pass alternately over and 
under successive wires. 

Pattern where each wire goes alternately over 
two wires and then under two successive wires. 

Wires running the long way as the cloth is 
woven. 



14. Weave 


The pattern of interlaced warp and shute wires 
determined by the sequence in which individual 
warp wires are raised and lowered by the 
heddle for passage of shuttle carrying and 
shuttle wire. 

15. Weft See 'shute'. 

16. Woof See 'shute'. 

NOMENCLATURE 

a = Surface Area to Unit Volume Ratio (1/0.3048 m) ( 1 /ft) 

B = Thickness (0. 3048 m) (ft) 

= Computed Pore Diameter (Inscribed Circle, (0. 3048 m) (ft) 

d = Wire Diameter (2. 54 cm) (in. ) 

d = Shute Wire Diameter (2. 54 cm) (in. ) 

s 

d w = Warp Wire Diameter (2. 54 cm) (in. ) 

N = Wires Per Inch (1/2. 54 cm) (l/in. ) 

N = Shute Wires Per Inch (1/2. 54 cm) (1/in. ) 

s 

N^ = Warp Wires Per Inch (1/2.54 cm) (l/in.) 

P = Percent Open Area 

Wioo = Weight, (kg/4. 22 m^) 100 ft^ Stainless Steel Mesh ( lbm) 
e = Void Fraction 

= Manufacturer's Absolute Micron Rating ( 1 fx = 0.0000394 in.) 

MESH SUPPLIER CODE 

1. Cambridge Wire Cloth Company 
Cambridge, Maryland (301) 228-3000 

2. Gerard Daniel and Company, Inc. 

New Rochelle, New York (914) 235-2525 

3. C. O, Jelliff Corporation 
Southport, Connecticut (203) 259-1615 
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4. TET/Kressilk 

Los Angeles, California (213) 283-3791 

5. Wire Cloth Enterprises, Inc. 

Pittsburgh, Pennsylvania (412) 731-1390 

6. Unique Wire Weaving Company, Inc. 
Hillside, New Jersey (201) 688-4600 


TWILLED DUTCH a 


N 

w 

N 

s 

b 

H-a 

Thickness (in. ) 

$/ft 2 

400 

2, 500 

5 

0. 0025 

300 

200 

1, 500 

13 

0. 0051 

50 

200 

1, 400 

14 

0. 0054 

41. 40 

200 

1, 200 

22 

0. 0052 

35. 90 

200 

1, 150 

12 

0. 0059 

34. 00 

200 

900 

20 

0. 0056 

28. 75 

200 

600 

28 

0. 0060 

22. 50 

150 

800 

21 

0. 0070 

26. 25 

120 

600 

34 

0. 0085 

18. 40 

120 

400 

45 

0. 0086 

10. 50 

120 

330 

50 

0. 0100 

10. 20 

120 

290 

59 

0. 0098 

8. 90 

120 

250 

60 

0. 0097 

8. 50 

120 

200 

70 

0. 0093 

6. 00 

120 

180 

70 

0. 0095 

5. 50 

120 

160 

80 

0. 0098 

5. 50 


cl- 

Products of Unique Wire Cloth Co., Inc. 
Based on bubble point test in alcohol 
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7 

• 50E-04 

2343 

.780 

1 .01E-03 

6.76E+03 

1 . 23E-03 

760/ .004 

4 

57.8 6.4 

6 

.67E-04 

2326 

.806 

1 . 06E-03 

5.55E+03 

1 .03E-03 


62# .0045 

4-, La 








52.0 

. 

CO 

7.50E-04 

2427 .773 

9.69E-04 

7 . 40E+03 

1 . 26E-03 

? 

64# .0045 

4 -,(* 
50.7 

8.7 

7.50E-04 

2511 .765 

9.27E-04 

8 • 09 E + 03 

1 . 28E-03 


65/ .0075 
765/ .0065 

! 1 3 

I I 

26.3 

'/T-.S 

33.4 

26.3 

19.3 

1 • 25E-03 
1 .08E-03 

2726 .574 
2660 .640 

6. 57E-04 
7.40E-04 

2.82E+04 
l • 87 E + 04 

3.79E-03 

2.65E-03 

7 

66/ .0045 

4 








49.4 

9.0 

7.50E-04 

2596 .757 

8.88E-04 

8 .82E+03 

1 .31 E-03 

766/ .0040 

4 








54.1 

7.1 

6 • 67E-04 

2573 .786 

9 .29E-04 

7.15E+03 

1 . 08E-03 


70/ .009 

13.7 

43.4 

1 .50E-03 

3119 .415 

4.40E-04 

8 . 46E+04 

8 . 7 0E-03 

770/ .008 

1 t, 3 . Co 

19.4 

33.2 

1 . 33E-03 

3025 .496 

5.24E-04 

4.96E+04 

5. 42E-03 

770/ .007 

1 ,7.3,G> 
26.0 

24.7 

1 . 17E-03 

2939 .571 

6.07E-04 

3 . 09 E+04 

3.57E-03 

770/ .006 

1 ,7,V 
33.6 

17.7 

1 .00E-03 

2862 .642 

6.90E-04 

1 • 99 E + 0 4 

2.42E-03 

770/ .0037 

54.9 

6*4 

6.17E-04 

2726 .790 

8.82E-04 

7.35E+03 

9 • 88E-04 
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P vV b a 

IOC 

42.0 9.7 6.67E-04 3517 

«-,G> 



47.9 

7.3 

5 .83E-04 

3471 

7 

90# .006 






21 .2 

23.8 

1 .00E-03 

3856 

790# .0055 

25.5 

19.7 

9.1 7E-04 

3786 

790# .005 

30.3 

16.0 

8.33E-04 

3721 

790# .0035 

4-> 

46.9 

7.5 

5 .83E-04 

3557 


94# .0040 4- 

36 

794# .0035 4-, 


38.9 10.4 6.67E-04 3786 

4> ,&> 

45.0 7.8 5.83E-04 3731 


100# .005 
7100# .0045 
7100# .004 
7100# .0035 
7100# .003 
7 

105# .003 


'P,3,U 

25.0 18.1 8.33E-04 4215 

30.3 14.4 7.50E-04 4134 

<A, S > 

36.0 11.1 6.67E-04 4060 

>,P*> 

42.3 8.4 5.83E-04 3994 

' i k 

49.0 6«l 5.00E-04 3936 


4 , 5 , <#> 

46.9 6.4 5.00E-04 4150 


120#. 004 l,X 

27.0 13*8 6.67E-04 5018 

7120#. 0037 

30.9 11 .6 6.17E-04 4950 

7 120#. 0026 

47.3 5.5 4.33E-04 4739 

7120#. 0035 2 

33.6 10*3 5.83E-04 4907 


130# .0038 I JL 

25.6 13.5 6.33E-04 5466 

7130#. 0034 

31.1 10*6 5.67E-04 5358 


620 3.58E-04 4.23E+04 1.47E-0 


A 

FT * 


M d. 

P w 

b 

a_ 

e 

D 

b O? 

b 



1 Oo 


c 

6^ 

<£p 

1 35# .0023 

4 -, 5 

47.5 

4.8 

3-83E-04 

5329 

.745 

4.26E-04 

1 .96E+04 

6.91E-S 4 

7 

1 40* .0033 

1 Z 









28.9 10.9 

5.50E-04 

581 4 

.600 

3.20E-04 

5. 16E+04 

1 . 53E-03 

?140> .0029 

1 7. / 3 < (° 
35.3 

8.2 

4.83E-04 

5696 

.6 56 

3.54E-04 

3.65E+04 

1 . 12E-03 

7 

1 45* .0022 

4.5, u 









46.4 

4.8 

3.67E-04 

5738 

.737 

3.91E-04 

2.22E+04 

6.75E-04 

7 

150* .0030 










30.3 

9.6 

5.00E-04 

6201 

.612 

3.06E-04 

5.13E+04 

1 . 33E-03 

? 150. .0026 









37.2 

7.0 

4.33E-04 

6070 

.671 

3.39E-04 

3.54E+04 

9.62E-04 

7 

160* .0028 

11 









30*. 5 

8.9 

4.67E! 04 

6610 

.614 

2.87E-04 

5. 40E+04 

1 • 24E-03 

? 1 60. .0025 

i 









36.0 

7.0 

4. 17E-04 

6497 

.662 

3. 13E-04 

4 . 02E+04 

9.52E-04 

7 

165* .0019 

A-,U, 

47.1 

4.0 

3.1 7E-04 

6519 

.7 42 

3.47E-04 

2.44E+04 

5.75E-04 

? 

170. .0026 

»,T- 









31.1 

8.1 

4.33E-04 

7007 

.620 

2.74E-04 

5.53E+04 

1 . 1 3E-03 

? 170. .0024 

35.0 

6.3 

4.00E-04 

6922 

.654 

2.90E-04 

4.43E+04 

9 . 35E-04 

7 

180* .0025 

• .~z- 

30.3 

8.0 

4. 1 7E-04 

7441 

.612 

2.55E-04 

6.1 5E+04 

1 . 1 IE-03 

7180. .0024 

~L 

32.3 

7.3 

4.00E-04 

7392 

.630 

2.63E-04 

5.50E+04 

1 .01E-03 

7180. .0023 

34.3 

6.7 

3.83E-04 

7344 

.6 48 

2.71E-04 

4.92E+04 

9 . 1 3E-04 

7 

190* .0022 

3 

33.9 

6*4 

3.67E-04 

7763 

• 6 44 

2.55E-04 

5.33E+04 

8 .84E-04 

7 

200* .0023 

1,1 

29.2 

7.5 

3.83E-04 

8299 

.602 

2.25E-04 

7 . 28 E+04 

1 i 06E-03 

7200. .0021 

1 








33.6 

6.2 

3 . 50E-04 

8178 

.6 42 

2. 42E-04 

5. 63E+04 

8 . 49E-04 

7200. .0016 

4,S / & 

46.2 

3.5 

2.67E-04 

7916 

.736 

2.83E-04 

3.08E+04 

4.92E-04 

7 

230* .0014 

46.0 

3.1 

2.33E-04 

9109 

.734 

2.46E-04 

3.59E+04 

4.33E-04 
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N d_ 


P VV b 

lOO 


a. 


-D. 


bo! 

e z 


FT* 


250# .0016 

S,£# 

36.0 

4.5 

2.67E-04 

10151 

.662 

2.00E-04 

6.28E+04 

6.09E-04 

? 









280# .0014 

S 

37.0 

3.8 

2.33E-04 

11338 

.669 

1 .81E-04 

6 • 70E+0 4 

5.21E-04 

?280j .0012 

S 

44. 1 

2.8 

2.00E-04 

11136 

.722 

1 .98E-04 

4.76E+04 

3.84E-04 

? 

325# .0012 

S 

37.2 

3.2 

2 . 00E-04 

13151 

.671 

1 .56E-04 

7.68E+04 

4.44E-04 

7325# .0011 

41 .3 

2.7 

1 .83E-04 

13012 

.702 

1 .65E-04 

6.30E+04 

3.72E-04 

? 

400# .0010 

5 

36.0 

2.8 

1 .67E-04 

16241 

.662 

1 .25E-04 

1 .00E+05 

3.81E-04 1 

?40Q# .0009 

5 

41 .0 

2.2 

1 .50E-04 

16027 

.699 

1 .33E-04 

7 .87E+04 

3.07E-04 
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"TWILLE.Q 3QUAR.E. 


M dl P 

\>y 

iOO 

to 

cl € 

r o 

0 

b oZ 

b 

A, 

FT* 

100* .0045 3 -Mr 
30-3 

13.7 

7 .50E-04 

3952 .630 

4.58E-04 

2.96E+04 

1 .89E-03 

£ Z.ZS 

?1 10* .0045 3,C» 
25.5 

15.2 

7.50E-04 

4387 .589 

3.83E-04 

4. 16E + 04 

2. 16E-03 

?1 20* .0037 3 

30.9 

11.1 

6 . 1 7E-04 

4737 .635 

3 • 86E-04 

3 • 43E + 04 

1 .53E-03 


7130* .0038 3 , U> 
25.6 

12.8 

6 .33E-04 

5184 .590 

3 -24E-04 

4 . 89E + 04 

1 . 82E -03 


? 1 40* .0033 3> 
28.9 

10.4 

5 .50E-04 

5546 .619 

3 • 20E -04 

4 . 42E+04 

1 .44E-03 


7150* .0028 3 

33.6 

7.9 

4.67E-04 

5894 .656 

3.22E-04 

3.77E+04 

1 .08E-03 


7160* .0028 3 

30.5 

8.5 

4.67E-04 

6321 .631 

2.87E-04 

4.68E+04 

1 • 1 7E-03 


165* .0026 3 

32.6 

y.5 

4 .33E-04 

6494 .648 

2.88E-04 

4 . 35E + 04 

1 .03E-03 


7170* .0026 

31 .1 

7.8 

4.33E-04 

6708 .637 

2-74E-04 

4.81E+04 

1 .07E-03 


7180* .0025 3, U 
30.3 

7.6 

4.17E-04 

7114 .630 

2.55E-04 

5.32E+04 

1 .05E-03 


7190* .0024 3 

29.6 

7.4 

4.00E-04 

7518 -624 

2.39E-04 

5 . 80E+04 

1 . 03E - 03 


7200* .0023 3,4- > 
29.2 

7.2 

3.83E-04 

7920 .621 

2.25E-04 

6.24E+04 

9 .96E-04 

$Z.Z5 

7200* .0020 5 

36.0 

5.4 

3.33E-04 

7830 .674 

2 .50E-04 

4 • 50E + 04 

7.34E-04 


7230*. 0018 3 

34.3 

5.0 

3 .00E-04 

9028 .661 

2* 12E-04 

5 . 59E+04 

6 . 86E - 04 


7250**0016 3,+ ,fc> 
36 .0 

4.3 

2.67E-04 

9788 .674 

2.00E-04 

5 • 63E+04 

5 .87E-04 


7270* *0016 3,4-, to 
32.3 

4.7 

2.67E-04 

10633 .646 

1 .75E-04 

7.24E+04 

6.40E-04 


7300* .0015 3. C* 
30.3 

A • 6 

2.50E-04 

11856 .630 

1 .53E-04 

8 • 87E+04 

6 .3 IE-04 

80 

7325* .0014 5,4-,£> 

29.7 4.3 

2.33E-04 

12857 .625 

1 .40E-04 

9 .87E+04 

5 .97E-04 


7400* .0010 3, 4- ( '5 , 
36.0 

2.7 

1 .67E-04 

15660 .674 

1 .25E-04 

9 .00E+04 

3 .67E-04 


7450* .00106 5 

27.4 

3.4 

1 .77E-04 

17880 .605 

9 . 69E-05 

1 • 54E + 05 

4 .82E-04 


7500* -0010 4 

25-0 

3.4 

1 .67E-04 

19962 .584 

8.33E-05 

1 .95E+05 

4 .88E-04 

$4-3.00 

7508* .00108 £ 

20.4 

4.1 

1 .80E-04 

20498 .539 

7.40E-05 

2.61E+05 

6.20E-04 

7508* .000866 S 
31 .4 

2.6 

1 .44E-04 

20037 .639 

9. 19E-05 

1 . 42E + 05 

3.54E-04 


7635*. 0008 4- 

24.2 

2.8 

1 .33E-04 

25395 .577 

6.46E-05 

2.58E+05 

4.01E-04 

$I3T00 

7850* .00063 4- 
21 .6 

2.3 

1 .05E-04 

34197 .551 

4.55E-05 

4 . 04E+05 

3.46E-04 

$ »3°I.OO 
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DUTCH 


NU Kl«. b 

W 5 W 5 100 

0- £ 

\.3 b 
£. z 

1.3 b $ 

Z 1 pt3 

/* c. 

10* 52* .028* .023 3 *• 



2.31E-02 $ 3.L0 


125.2 6.17E-03 
712*64* .023* .0165 3 

827 .589 

1 .58E+04 

31 S 

82.9 4.67E-03 
712*88* .014* .013 3,Z 

973 .641 

1 .40E+04 

1 .48E-02 

S3S 

58.2 3.33E-03 
714*64. . 020* .0165 3 

1292 .647 

1 . 73E+04 

1 .04E-02 

330 

81.4 4.42E-03 

71 4*88* .020* .013 3 

1043 .627 

1 .59E+04 

1 .46E-02 

320 

71.1 3.83E-03 

714*95* .015* .012 3 

1262 .625 

2.03E+04 

1.28E-02^S.T5 

30S 

57.5 3.25E-03 

714* 100* .015* .012 3 

1381 .642 

1 .95E+04 

1 . 02E-02 


60.0 3.25E-03 
720* 120* .014* .010 3 

1443 .627 

2.24E+04 

1 .08E-02 

ZSo 

55.7 2.83E-03 
724* 1 10* .015* .0105 1 ,3 ,4-, 

1786 .602 

3 • 24E+04 

1 . 02E-02 
1.19E-02 

zoo 

63.5 3.00E-03 
724* 175* .01 1* .0065 3 

1795 .572 

3.84E+04 


35.3 2.00E-03 
730* 1 50* .009* .007 ~l ,4- 

2353 .643 

3.49E+04 

6.30E-03 

140 

34.2 1.92E-03 
730* 160, ,009* .007 Z.,3, C» 

2352 .639 

3 • 38E+04 

6. 1 IE-03 

\ZO 

36.0 1.92E-03 

740*200* .007* .0055 

2479 .620 

3.98E+04 

6.48E-03 

28.2 1 . 50E-03 

750*250* .0055* .0045 

3162 .619 

5.08E+04 

5.08E-03 

13 

23.2 1.21E-03 

•> 

* SUPPLIER. CODE! 

3985 .611 

6.67E+04 

4.20E-034 Z.4Z. 
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REVERSE. PL&IW DUTCH 


W ,oo ^ 

a. 

£ 

13 b a? 

1.3 b $ 

£ Z FT. 2, 


1 5/ 64* .0276/ .01 59 Z* 

99.5 4.95E-03 

1034 .593 

1 .95E+04 

1 .83E-02 


718/64/ *0276* .0159 Z. 

110.8 4.95E-03 

1136 .547 

2.77E+04 

2.15E-02 


? 18. 107* .0238/ .00943 Z 

68.8 3.56E-03 

1418 .608 

2.51E+04 

1 .2SE-02 


720/ 107/ .0238* .00943 Z 

73.6 3.56E-03 

1 49 1 .581 

•3.04E+04 

1 .37E-02 


7 36/ 173* .0138/ .00593 l^Z 

46.3 2.14E-03 

2579 .562 

5.86E+04 

8.8 IE-03 


738/ 132/ .0159# .00793 \,Z 

67.3 2 • 65E-03 

2393 .486 

8 . 36E+04 

1 • 46E-02 


7 41/ 173# .01 19/ .00593 \,Z 

43.2 1.98E-03 

2790 .559 

6.41E+04 

8.24E-03 


7 48/ 173/ .01 19/ .00593 \,Z 

47.9 1.98E-03 

3049 .511 

9. 18E + 04 

9.87E-03 


7 50/ 17 5/ .012/ .006 4- 

50.7 2.00E-03 

3161 .487 

1 . 10E+05 

1 . 10E-02 $ Co. 7.0 


? 50/ 180/ .012/ .0057 4- 

48.5 1.95E-03 

3174 .496 

1 . 04E+05 

1 • 03E-02 

100 

769/292/ .0087/ .00343 1 ,Z 

31.4 1.30E-03 

4619 .510 

1 .38E+05 

6 . 47E-03 


770/280/ .0083# .0035 4- 

30.0 1.27E-03 

4564 .524 

1 • 26E+05 

6.03E-03 $ 3. 'SO 

S4- 

776/292/ .0079# .00343 1,7. 

30.0 1.23E-03 

4909 .506 

1 . 50E+05 

6-24E-03 


? 102*635* .0064* .001 58 l,Z 

20.1 7.97E-04 

7676 .489 

2.55E+05 

4.33E-03 


7114/635/ .0051* .00158 1 ,Z 

16.0 6.88E-04 

8409 .529 

2.26E+05 

3.20E-03 


? 120* 180. .004/ .0035 Z 

15.8 9.17E-04 

4550 .651 

5.82E+04 

2.81E-03 


? 120*200* .004. .0031 7. 

14.4 8 . 50E-04 

4785 .658 

5.85E+04 

2.55E-03 


7120/400/ .0036* .0022 7. 

12.7 6.67E-04 

7088 .616 

1 . 15E+05 

2.29E-03 £ 4.Z0 


7 120/670/ .0047/ .0015 4- 

14.7 6.42E-04 

8894 .538 

2.28E+05 

2.38E-03 ^ ! ! .4-0 

LI' 

? 125/600/ .0049/ .0016 4- 

16.1 6.75E-04 

8608 .517 

2.43E+05 

3.28E-03 i U. CO 

Z1 

7125/ 400/ .0036/ .0022 Z 

13.0 6 . 67E-04 

7243 .606 

1 .24E+05 

2 • 36E-03 


? 125/600/ .0048/ .0016 Z. 

15.7 6.67E-04 

8621 

.524 

2.34E+05 

3.1 5E-03 


7130/635/ .0051/ .001 58 \7. 

17.7 6.88E-04 

9091 .479 

3.22E+05 

3 • 90E-03 


7 1 40/720/ .0043/ .001 4 4- 

14.2 5.92E-04 

10044 .514 

2.94E+05 

2.91E-03 £ il.00 

Z\ 

71 35/635/ .0051 / .001 58 1 ,7_ 

18*3 6.88E-04 

9308 .463 

3.61E+05 

4. 17E-03 


71 52/747/ .0044/ .00134 1 % Z 

15.3 5.90E-04 

10635 .476 

3.83E+05 

3.38E-03 


71 55/850/ .004/ .0012 4- 

13*1 5.33E-04 

11369 .503 

3. 55E+05 

2.74E-03 $ I4-.40 

7.0 


o 
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TWILLED OtlTXH 


1,3 b cl L3b % 

N s d w a s ^00 k ^ 6 ^ & 1 PtT h-o- 


16*200* .01 1 * *010 3 -* 

72.7 

2.58E-03 

2713 .430 

1 .33E + 05 

1 .81E-02 

130 

720*250* *010* *0085 2> <*-,<!= 

66.8 

2.25E-03 

3345 .399 

2.06E+05 

1 .84E-02 

80 

720*250* *0036* .008 Z 

53.3 

1 .63E-03 

4036 .340 

2.99E+05 

1 .84E-02 


720*280* .01 1* .0075 3 

60.5 

2.17E-03 

3465 .435 

1 • 78E + 05 

1 .49E-02 

132. 

720*200* .01 1* .010 3 

75.3 

2.58E-03 

2803 .411 

1 • 5 7E+05 

1 .99E-02 

no 

720*350* .007* .006 3 

44.6 

1 .58E-03 

4514 .430 

2.27E+05 

1 .1 IE-02 

Ld8 

728*500* .007* .0045 3,<o 

38-0 

1 -33E-03 

5897 .423 

3.37E+05 

9.68E-03 

CoO 

730*250* .010* .008 Z. 

65.0 

2-17E-03 

3533 .393 

2.28E+05 

1 .82E-02 


730*260* .010* .008 3> 

67.2 

2.17E-03 

3657 .372 

2 • 72E+05 

2.03E-02 

25 

730*370* .010* .006 ^ 

55.1 

1 -83E-03 

4526 .391 

3 • 1 9E+05 

1 • 56E-02 

IOO 

730*500* .008* .006 3 

66.9 

1 • 67E-03 

6347 .188 

2.47E+06 

6.14E-02 

51 

740*570* .0071* .0040 4- 

37.4 

1 .26E-03 

6669 .398 

4.58E+05 

1 .03E-02 

lo 

750*700* .006* .003 3 

27.2 

1 .00E-03 

7858 .450 

3.97E+05 

6.42E-03 

SO 

780*700* .004* .003 1^3 

26.0 8.33E-04 

9702 .369 

7.51E+05 

7.98E-03 

^1,50 43 

7150*800* .0027* .0021 3 

16.3 

5.75E-04 

12458 .427 

6 • 35E+05 

4.09E-03 

35 

7164*800* .0028* .002 l 

15.9 

5 .67E-04 

12626 .432 

6 .31E+05 

3.96E-03 

Z.L> 

7165*800* .0028* .0020 

16.0 

5 .67E-04 

12654 .430 

6.38E+05 

3.98E-03 

2.1® 

7165*800* .0029* .0020 Z 

16.3 

5.75E-04 

12606 .426 

6.55E+05 

4. 12E-03 


7165* 1200* .0028* .0016 Z 

15.3 

5.00E-04 

16389 .381 

1 .20E + 06 

4.47E-03 


7165* 1400* .0028* .0016 i ,7. ( 4- 

17.1 5.00E-04 

18636 .306 
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APPENDIX B 


FLOW TEST DATA CORRELATION 
Screen Flow-Through Loss Data 

The screen flow-through loss data, plotted in Armour and Cannon form, 
are shown for all 10 screens in Figures 36 through 42. Also shown in these 
figures are data from all other known sources. These data have been 
normalized to agree with the geometric characteristics of our screens 
(especially specified pore diameter) as shown in Table VIII. 

Figure 36 shows the available data for the 325 x 2,300 screen, including 
gas data from the MDAC-MSFC contract NAS8-27685 (ref. 12) which agree 
with our data rather well. Water test data from Kressilk (ref. 13) and 
Wintec Corp. (ref 14) are also shown, together with the generalized Armour 
and Cannon correlation for all screens, which overpredicts the friction 
factor and thus the pressure loss for this screen by 170%. , 

Many data have been accumulated for the 200 x 1,400 screen, and our 
data are generally centered on the normalized available data, as shown in 
Figure 37. Some of the data scatter evident in Figure 37 may be caused by 
unaccounted differences in the screens tested (wire diameter, void fraction, 
etc. ) or to the rather large variations in the properties of the test fluids. 
Again, the generalized Armour and Cannon correlation is substantially above 
the correlation for our data. 

Figure 38 shows the correlation for the 720 x 140 reverse Dutch screen, 
which lies above, but close to the generalized Armour and Cannon correla- 
tion. Previous MDAC data also agree well with our data. Figure 39 shows 
the correlation for the 16 5 x 800 screen. The GDA data were definitely 
obtained with a different type of 1 6 5 x 800 screen (different wire diameter), 
but when normalized to our 165 x 800 screen characteristics, agree very 
well with our data. Figure 40 shows the correlation for our 50 x 250 plain 
Dutch screen, which again agrees very well with the normalized GDA data 
(ref. 15) using GH 2 and GN 2 . Figure 41 shows the correlation for the 
24 x 110 plain Dutch screen, which coincidentally falls right on the general- 
ized Armour and Cannon correlation, as does the Kressilk data. Therefore, 
the a and (3 values from the Armour and Cannon correlation were retained 
for the 24 x 110 screen. 

The correlations for the square weave screens are all shown in 
Figure 42, and were discussed previously. 

Channel Flow Loss Data 

2 ' 

The data, plotted as head loss in meters of T4. 5 N/cm LH 2 versus chan 
nel flow velocity in meters /sec with length (L) over channel height (s) as a 
parameter, are shown for the 10 screens in Figures 43 through 52. As 
anticipated in test planning, the pressure drop was linear with length: in 
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Figure 40. Flow Loss Correlation — 50 x 250 








Figure 42 . Flow Loss Correlation — Square Weave Screens 
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Figure 43. Channel Flow Pressure Drop Correlation — 325 X 2300 
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PRESSURE DROP - METERS OF 34.5 N/CM c LH, 
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Figure 45. Channel Flow Pressure Drop Correlation — 500 X 500 
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Figure 46, Channel Flow Pressure Drop Correlation — 720 X 140 
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Figure 47. Channel Flow Pressure Drop Correlation — 165 X 800 
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PRESSURE DROP - METERS OF 34.5 N/CM c LH 
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Figure 50 . Channel Flow Pressure Drop Correlation — 24 X 110 
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PRESSURE DROP - METERS OF 34.5 N/CM‘ LH, 







126 


DROP - METERS 0 







Figures 43 through 52, the plain symbols are for the 33 cm (13-inch) length, 
the primed symbols are for 66 cm (26-inch), and the double-primed symbols 
for the 99 cm(39-inch) length (except for the 325 x 2, 300 and 500 x 500 
screens where they correspond to 99, 132, and 198 cm (39, 52, and 
78-inches) respectively. Because the channel flow lies in the turbulent 
regime, the pressure drop varies essentially as the channel flow velocity 
squared, as shown by the slope of the lines on Figures 43 through 52. In 
general, the pressure drop increases with increasing L/ s , but in some cases, 
extreme channel height values cause jumps in the L/ s correlation because 
of the effect of channel height on friction factor; see, for example Figures 44, 
46 , 50 , and 52. 

To determine the effects of channel height (or effectively, e/D^) on our 
channel flow loss and the appropriate roughness parameter, our data were 
plotted on the Moody graph, as shown in Figures 53 through 62. In these 
figures, the collection of data points at R = 10^ or higher were the 
34.5 N/cm^ (50 psia) LH 2 data, while the single points at R = 5,000 were 
taken with ambient GN 2 . In our experiments, the direction of fluid flow was 
always in the direction of the shute wires, which should minimize the 
pressure drop due to the construction of the Dutch- weave screens. It was 
noted from the data that the f for the Dutch-weave screens was about half 
that for square- weave screens of about the same wire diameter and essen- 
tially identical geometry and flow conditions. This is shown, for example, 
by comparing the 24 x 110 screen (Figure 60) with the 40 x 40 screen 
(Figure 62), or the 1 6 5 x 800 screen (Figure 57) with the 150 x 150 screen 
(Figure 59). 
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Figure 53 . Dimensionless Channel Flow Loss Correlation — 325 x 2300 
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Figure 55. Dimensionless Channel Flow Loss Correlation — 500 x 500 
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Figure 56. Dimensionless Channel Flow Loss Correlation — 720 x 140 
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Figure 57. Dimensionless Channel Flow Loss Correlation — 165 x 800 
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Figure58 . Dimensionless Channel Flow Loss Correlation — 50 x 250 
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Figure 59. Dimensionless Channel Flow Loss Correlation — 150 x 150 


RELATIVE ROUGHNE 



CR 140 


VALUES OF (VD") FOR WATER AT 60*F (VELOCITY IN tdT X DIAMETER !N INCHES) IQCOG 

0.2 04 OS 08 I 2 4 6 0 10 20 43 60 60 100 200 400 600 801 OOP 2000 4-00 0 60 03 60031 

T ; f~~ f VALUES CF (VD’i'fDR ATMOSPHERIC AIR AT 60'E 1 T~f edcxTl ” i I 8QOOO I I 

4 6 8 10 20 I 40 [ 60 MOO. I 200 ] 400 1600 POO 003 2033 4COO 6000 | 10,000 2Q000j 4-3.000 [60,000 100,000 

- i — i 1 "i T |-j j f’-pi r r f — S — f- -r* r 1 -t v ri ll T . ' 7 "T t' — ‘ t — i 1 ""i J r rfi'i"rh-n L — n — i’-j i‘ i i' i i i i m r • i — r*-r‘ i‘ r I'rn-hT 


nf:i k s 1 n t \vrr\-\ nn- 




I i i ! ! r 1 ! 


!^n'fc^]TR/^srn<y^x_^ l 11 Li 1 ’ i I 7LP i ’ J !. Li JI'J. i. I Mill 

-jFUJ>V-^2Q\E-^ZCn^E — 1 — i — COMPLETE 1 TU R BULENCE, ; EOUoH; PIPES -f 


| Li J l J. i ! J J u |. i j j j ' 


a- i j irk 

lp~-l' Prtp:: 

ltvtz.' ~ Z i Ti ria<.'iTy 
M-poV—t — 


.. .. ..... . ■ j - | ; rfj-- 


\ ; I : ! ' i 


*\ 

......... . V 




i '~ r j iTi ii~i : iii 

1 ; j | ! . i! M ! i . - 04 

- i Lt---= r .03 


. 00 ! 

.0003 

.0006 


24 x no 


II I I I 1 ! II I I M I I : I ; Mill | I I I 1.. ; I - ■ ; I i ; ; I L_i..L j_LL • i.LITSSw 

^ ECiOC 3 4 5 6 6 |q4 2(10“) 3 4 5 6 0 | Q 5 2(10-’,' 3 4 5 6 8 ;q6 20O 4 ) 3 4 5 6 8 10 7 

REYNOLDS NUMBER ( V iN-g , D in ft, V iN^f ) 


-4^: — ~r-T> .000,0! 

5jOf^3^4 5 6 3 .qS 

5 '* °°Q Oo/^"' °°0, o 0 


Figure 60. Dimensionless Channel Flow Loss Correlation — 24 x 1 10 
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Figure 61 . Dimensionless Channel Flow Loss Correlation — 60 x 60 
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Figure 62. Dimensionless Channel Flow Loss Correlation — 40 x 40 




APPENDIX C 


EVALUATION OF SIDE-WALL CONTRIBUTION TO CHANNEL 

FLOW PRESSURE LOSS 

In the channel flow apparatus, the direction of flow is in the direction of 
the shute wires in the screen. It is assumed that the equivalent roughness of 
the screen is the ratio of the shute wire diameter to the hydraulic diameter 
for square weave screens, and the ratio of half the shute wire diameter to the 
hydraulic diameter for Dutch-weave screens. Because the other channel 
walls are made from rolled stainless -steel sheet, it is further assumed that 
the equivalent roughness of the other channel walls is approximated by the 
roughness of drawn tubing or e = 0. 000152 cm (0. 00006 in. ). Finally, it is 

assumed that the total roughness of the channel is the sum of the roughness 
contribution of each of the channel walls, or: 


e 

D h 


TOTAL 


<2 

- D h 


SCREEN \ 2w + 2s 


( 2w + 2s ) ' + ° h 


SIDEWALLS 


( 2s\ 

\ 2w + 2s / 


'h 


WALL 


/ IX \ 

\2w + 2s / 


(C-l) 


The sidewall contribution to the total is 



e 

1 

/ 2s \ 



D, 1 
h 

SIDEWALLS ' 2w + Zs ' 


e 

1 / w \ , e 

1 1 

/ 2s \ , e 

/ w \ 

D h 

SCREEN V 2w+2s / D h 

SIDEWALLS \ 2w + Zs ' D h 

WALL V 2w + 2s / 


or 


e SCREEN 


0. 000152 (2s) 

w + 0. 000152 (2s) + 0. 000152 w 


(C-3) 


The sidewall contribution from equation (C-3), based on the screens 
tested in the as -built channel height (s) and width (w) is shown in Table XIII. 
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TABLE XIII. - SIDEWALL CONTRIBUTION TO CHANNEL FLOW LOSS 


Screen 

Roughness 

e(cm) 

Actual channel 
spacing depth 

Sidewall 

contribution 

(%) 



(cm) 

(in. ) 

40 x 40 

0. 0254 

0. 455 

0. 179 

0. 043 



0. 554 

0. 218 

0. 053 



1. 034 

0. 407 

0. 098 



1. 133 

0. 446 

0. 108 

24 x 110 

0.01334 

0. 457 

0. 180 

0. 082 



0. 531 

0. 209 

0. 096 



1. Oil 

0. 398 

0. 182 



1. 085 

0. 427 

0. 195 

720 x 140 

0. 00546 

0. 622 

0. 245 

0. 269 



0. 747 

0. 294 

0. 323 



1. 514 

0. 596 

0. 652 



1. 638 

0. 645 

0. 705 

200 x 1,400 

0. 00203 

0. 635 

0. 250 

0. 703 




0. 300 

0. 842 



1. 499 

0. 590 

1. 643 



1. 626 

0. 640 

1. 780 

o 

ON 

o 

0. 01905 

0. 587 

0. 231 

0. 074 



1. 252 

0. 493 

0. 158 



1. 948 

0. 767 

0. 246 



2. 614 

1.029 

0. 330 

50 X 250 

0. 00572 

0. 589 

0. 232 

0. 243 



1. 234 

0. 486 

0. 508 



1. 963 

0. 773 

0. 799 



2. 609 

1. 027 

1. 069 

150 x 150 

0. 0066 

0. 612 

0. 241 

0. 220 



1. 278 

0. 503 

0. 458 



1. 974 

0. 777 

0. 706 



2. 639 

1. 039 

0. 942 

165 x 800 

0. 00254 

0. 610 

0. 240 

0. 548 



1. 255 

0. 494 

1. 121 



1. 984 

0. 781 

1. 761 



2. 629 

1. 035 

2. 320 

500 x 500 

0. 00254 

1. 090 

0. 429 

0. 975 



1. 643 

0. 647 

1. 463 

325 x 2,300 

0. 00127 

2. 890 

1. 138 

4. 71 



3. 444 

1. 356 

5. 563 
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APPENDIX D 


VISCOUS PRESSURE DROP IN A SPHERICAL ANNULUS 


The viscous pressure drop for steady incompressible flow between a 
spherical screen annulus was determined based on the approach given, for 
example, by Bird, Stewart, and Lightfoot, (ref. 27). Creeping flow was 
assumed so that inertial terms could be neglected in the momentum equation. 
It was also assumed that v r = v<& = 0 and vq =-. vq (r, 0). The coordinates and 
geometry are shown in Figure 63. 

The governing equations are: 

Continuity 


-1+ A. 


~ + — ?+ r— ( p r v ) 
St 2 5r r 
r 


1 


r sin 


0 60 (pv 0 


sin 


0 ) 


r sin 0 6 


V 


0 (D-l) 


CR 140 



Q 


Figure 63- Laminar Flow in a Spherical Annulus — Coordinates and Geometry 
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Momentum (9 - component) 
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' r* 


2 . 2 . 
r sin 9 


2 cos 9 

2 . 2 . 
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6V (6\. 

TT/ pg e 


Where 


1 8/2 6 \ 

“ “P( r TT/ 4 " 


r sin 9 69 


7( si 


sin0 7) + 


2 2 
r sin 9^5^ 


(£) 


The continuity equation simplifies to 


89 (v 9 sin 9) 


0 


Thus, 

Vq sin 9 = u (r) 


The momentum equation simplifies to: 


r 89 


1 6 ( 2 6v e\ 1 6 / 6 M 

T5r l r TP 4 2“ 7 68 V‘ m0 Tr) 


r sin 9 


6 l ■ a 6v 0 

60 I s m 0 TF 


V 


9 


2 72 - + pg e 

r sm 9 


Since 


5v r 


sin 9 


6 9 


X 0 

4" (v 9/ in 9) -v Q cos 9 


- Vq COS 9, 


(D-2) 


(D-3) 


(D -4) 


(D-5) 


(D-6) 


(D -7) 
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the terms given below can be shown to be zero; i. e. 


6 ( 6v 0 

sin 0 _~V sin9 —j -v Q = 


sin Q__(_ Vq cos 9) - v g 


2 6v 9 

= Vq sin 9 - sin 9 cos 9 -^q— - Vq 


2 6v 9 

Vq ( 1 - cos 9) - sin 9 cosO^q- - Vq 


D-8) 




cos 9 ( cos 9 Vq - sin 9 -^q- 


cos 9~ftQ~( v Q s i- n 9) = 0 


With sin 0 = u (r). 


A _ 1 6 P /I .Id / 2 du\\ 

r &9 ^ ysin 9 ^2 dr \ dr / J 


(D-9) 


Using separation of variables gives 


sin 9 |P = B 


;d-io) 




:d-id 


Integrating equation (D-10) over the limits Pp P^, and 9p 9 gives 

P Q 


dp = 


9, 


Bd9 
sin 9 


(D - 1 2) 


AP = P - P 2 = - B In 


tan 9^/2 


tan 9^/2 


( D- 1 3 ) 
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Integrating equation (D-ll) twice gives 



M 2 dr)= > r+C 

(D- 14) 


2 du Br^ . 

r -j — - + C 

dr 2 |jl 

(D-15) 

and 

Br c * n 

2 jjl r 

(D-16) 

The boundary conditions 

are 



u = 0 at r = R 

u = 0 at r = KR 

( D- 17) 

Thus, the integration constants in (D-16) are obtained from 



0 = + D 

2g. R 

(D-18) 

and 

_ BKR C „ 

0 _ - + D 

2 jjl KR 

( D-19) 

Solving the above gives 

B (R - KR) 

(D-20) 


< k + i > 

(D-21) 

Thus, u (r) is 

“ - -f t -*♦«(> 4) 

(D-22) 
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or. 


u = 


APR 

1 - 

1 * K (‘ 

-f). 

2 p 

In 

"tan 92/21 


.tan 9 / 2 J 


-Vg sin 9 


The total flowrate is, at any point, 

m = f pv Q dA cg or Q ~ J Vg dA 

Thus, at 0 = tt/2, dA = 2fTrdr and 

u = sin 9 = v 

therefore, 

• f R 
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The flowrate is 
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The pressure drop is 


or, since 


AP 


6 p Q 

ttR 3 (1 - K) 3 



s = (R - KR) 


AP 


6 p Q 


IT S 


In 


tan 0^/2 
tan 9^/2 


(D -28) 
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APPENDIX E 


WALL-SCREEN SPACING AND TANK OPTIMIZATION ANALYSIS 
Wall-Screen Spacing Analysis 

Each tank was analyzed with all 10 screens and with annulus gaps from 
1% to 5% of tank volume. The results of this analysis for the 5,000/4 tank 
are shown in Figures 64 and 65. Figure 64 shows the performance of the 
10 screens in terms of safety factor and puddle residual for the 1% annulus 
gap (the specified minimum). Since the avowed purpose of the ultimate 
tradeoff study was to minimize system weight, it was logical to select screens 
with adequate performance at the minimum annulus gap. For system/screen 
comparison, a safety factor of 2 was chosen as representative of adequate 
performance. Figure 64 shows that at a 1% annulus gap, only the 325 x 2, 300 
and 200 x 1,400 screens gave a safety factor of 2. Since the 325 x 2, 300 
screen was substantially lighter than the 200 x 1, 400 screen and had a 
higher bubble point, it was the logical choice. 

Figure 65 shows that at 2% annulus gap, the 150 x 150 screen had com- 
parable performance to the 325 x 2, 300 screen in the 1% annulus. 

Figure 66 shows the results for the 500/4 tank. At a 1% annulus gap, 
five screens had adequate performance, of which the 325 x 2, 300 screen was 
the lightest in weight (except for the 500 x 500 screen). To use the 150 x 150 
screen, with a safety factor of 2, the annulus gap would have had to be 
increased to 1.27%, which would have incurred an increased residual of 
2. 7 kg (6 lb), while saving screen weight of 6. 8 kg (15 lb), for a net weight 
reduction of 4. 1 kg (9 lb). Similarly, Figure 67 indicates that eight screens, 
including both the 325 x 2, 300 and 150 x 150, would have adequate perform- 
ance in the 500/2 tank, while Figures 68, 69, and 70 indicate that all 10 
screens would have adequate performance for the remaining tanks. The 
weight savings using the 150 x 150 screen instead of the 325 x 2, 300 screen 
ranged from 5. 9 kg to 1. 1 kg (12. 9 lb to 2. 5 lb) (which were not significant 
weight savings from a system standpoint) while sacrificing an order of 
magnitude in bubble point. 

The analysis described above was for the maximum outflow rate of l%of 
tank volume /minute. At 0.01% outflow rate (or at 0. 1% TVS flow rate) there 
was absolutely no sensitivity to screen type or annulus gap; all 10 screens had 
adequate performance at the minimum 1% annulus gap. It is clear that at 
this g-level, flowrate, and annulus gap, the tankage systems flow character- 
istics were such that a meaningful tradeoff analysis in terms of residual, 
annulus gap, and performance was not possible. 

Therefore the analysis was extended to examine the effects of high out- 
flow rate (3% tank volume /minute) and high g-levels (10"^ to 10“^ g's) on 
system performance. Only the high L/D tanks (5,000/4, 500/4, and 50/2) 
were analyzed because these were generally the more severe cases. 

Figure 71 shows the effect of increasing the flowrate to 3% /minute on the 
performance of the 325 x 2, 300 screen in the 5,000/4 tank. The annulus gap 
had to be increased to 2% to achieve adequate performance, and even then, 
increased residual resulted. In Figure 72, the performance effects for the 
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Figure 65- Screen Performance in 2% Annulus for 5,000/4 Tank 
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Figure 67. Screen Performance in 1% Annulus for 500/2 Tank 
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Figure 70. Screen Performance in 1% Annulus for 50/1 Tank 
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Figure 72. High Outflow Rate Performance of 150 x 150 Screen in the 5,000/4 Tank 
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150 x 150 screen in the 5, 000/4 tank are shown. The annulus gap had to be 
increased from 2% to 5% to obtain adequate performance. Figure 73 shows 
that for the 500/4 tank with increased flowrate, the annulus gap had to be 
increased from 1% to 2% for the 325 x 2, 300 screen, and from 2% to 3% for 
the 150 x 150 screen. (These were the only screens analyzed further since 
they represented the minimum weight and maximum bubble-point screens. ) 
Figure 74 shows that for the 50/2 tank, the 325 x 2, 300 screen had adequate 
performance with a 1% annulus gap at increased flowrate but that the annulus 
gap for the 150 x 150 screen had to be increased from 1% to 2%. 


A similar trend occurred when the g -level was increased from 10“^ g's 
to 10”^ g's. Figure 75 shows the performance of the 325 x 2, 300 screen in a 
1% annulus in the 5, 000/4 tank at various g-levels. At 10“^ g's, the 
325 x 2, 300 screen no longer had adequate performance, and the annulus gap 
had to be increased. In Figure 76, the annulus gap is 2% in the 5,000/4 tank. 
The 325 x 2, 300 screen had adequate performance at all g-levels, but the 
150 x 150 screen did not have adequate performance at even 10”^ g's. If the 
annulus gap was increased to 3%, however, the 150 x 150 screen would have 
adequate performance at even 10“^ g's. At 10“^ g's, the hydrostatic head in 
the 5, 000/4 tank exceeded the bubble point of the 150 x 150 screen by a factor 
of 3, so that it could not practically be used. Similarly, in Figure 77 the 
performance of the 325 x 2, 300 screen in a 1% annulus and the 150 x 150 
screen in a 2% annulus is shown for the 500/4 tank, and in Figure 78, the 
performance for these screens in the 50/2 tank is shown. 


Pump /Standpipe Optimization 

The standpipe optimization equation (equation 44) is plotted as the nearly 
horizontal lines in Figures 79 through 90 for the six tanks and two values of 
TVS flowrate. 


From equations (42), (45), and (55), the overall efficiency as a function 
of the total fluid power, including the standpipe loss, was determined and 
plotted in Figures 79 through 90. Where these lines crossed the standpipe 
optimization lines, the intersection was the value of standpipe diameter which 
gave the minimum system weight for that annulus gap, screen, and standpipe 
material. It was assumed that either stainless steel or aluminum could be 
used for the standpipe because of the presence of the slipjoint in the stand- 
pipe (Figure 1) which allowed differential thermal expansion. With the design 
value of standpipe diameter as an input, and with the other head losses known, 
the total system weight analysis was completed. 

System Weight Optimization Analysis 


The overall system weight was divided into four categories: annulus 
residual, standpipe residual, pump boiloff, and hardware, (which includes 
standpipe weight, screen weight, pump and motor weight. ) The four cate- 
gories above plus the total weight is shown for the six tanks, for 30 and 
300-day missions, for various screens and standpipe materials versus 
annulus gap in Figures 91 through 103. The results for the 30 -day mission 
are shown in Figures 91 through 96 and indicated that minimum system 
weight was achieved with the minimum annulus gap because of the strong 
influence of annulus residual. 
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Figure 77. High G-Level Performance of Both Screens in the 500/4 Tank 
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Figure 791 Standpipe Optimization at 1%/Min TVS Flow in the 5,000/4 Tank 
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Figure 80. Standpipe Optimization at 1%/Min TVS Flow in the 500/4 Tank 
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Figure 81. Standpipe Optimization at 1%/Min TVS Flow in the 500/2 Tank 
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Figure 83. Standpipe Optimization at 1%/Min TVS Flow in the 50/2 Tank 
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Figure 84. Standpipe Optimization at 1%/Min TVS Flow in the 50/1 Tank 
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Figure 85 ! . Standpipe Optimization at 0.1%/M in TVS Flow in the 5,000/ 4 T ank 
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Figure 86. Standpipe Optimization at 0.1%/Min TVS Flow in the 500/4 Tank 
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Figure 87 . Standpipe Optimization at 0.1%/Min TVS Flow in the 500/2 Tank 
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Figure 88. Standpipe Optimization at 0.1%/Min TVS Flow in the 500/1 Tank 
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Figure 89- Standpipe Optimization at 0.1%/Min TVS Flow in the 50/2 Tank 
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Figure 91. Weight Optimization for 5,000/4 Tank for 30-Day Mission 
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Figure 93> Weight Optimization for 500/2 Tank for 30-Day Mission 
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Figure 97. Weight Optimization for 5,000/4 Tank for 300-Day Mission 
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Figure 98. Weight Optimization for 500/4Tank for 300-Day Mission 
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However, for the 300 -day mission, the pump boiloff became very 
important, as shown in Figures 97 through 102, so that an optimum annulus 
gap was found, at which minimum weight occurs. For the 5,000/4 and 
500/4 tanks, this optimum was at about 2.0 to 2. 2% annulus gap for both the 
325.x 2, 300 and the 150 x 150 screens. For the 500/2 and 50/2 tanks, the 
optimum gap was at about 1. 6%, and for the 500/1 and 50/1 tanks the opti- 
mum was at an annulus gap of about 1.4 to 1. 5%. In Figure 103, for the 
5, 000/4 tank at 0. 1% tank volume/minute TVS flow, the optimum again 
occurred at the minimum annulus gap, as was the case for the other five 
tanks. 
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APPENDIX F 


PROPELLANT HEATING FROM ELECTRIC PUMPS 


Most of the power input to the electric motor/pump in the LH? tank is 
lost in electric motor inefficiency, principally windage losses ana friction, 
which directly heat the LH 2 . Some of the remaining power to the pump is 
lost as pump inefficiency, principally friction losses, with the remaining 
power being used as fluid power. However, nearly all of the fluid power is 
used to overcome flowing friction and only a relatively small fraction is 
stored in the LH 2 as fluid kinetic energy. All other power loss is essentially 
dissipated in the LH 2 as heat, leading to "boiloff. " 

To determine the maximum energy storage capability of the system, it 
can be assumed that over a long period of time the pump imparts momentum 
to not only the fluid in the annulus, but to all of the fluid in the tank by 
momentum exchange. The assumption that all of the fluid in the tank 
acquires the maximum velocity of the fluid in the annulus at the specified 
TVS flowrate gives the maximum energy storage capacity of the system. 

That energy is, for the 141. 6 m^ Q, 000 ft^) tank, simply the mass of liquid 
times the maximum dynamic head in the annulus at a TVS flowrate of 1% 
of tank volume /minute, or: 


K. E. 


V P H d 


141.6 m 


3 


(64. 08 kg/m 3 ) 0.0488 m 9.807 


watt- sec 
kg - m 


= 4, 335 watt- sec 

Compare this to the total energy entering the tank through the pump in 
30 days: 

9.7 5 watts x 2, 592, 000 sec = 25, 250, 000 watt- sec. 

The ratio 4, 335/25, 250, 000 = 0. 017% of energy input remains stored as 
kinetic energy for a 30-day mission, or 99. 983% is dissipated as heat to the 
LH 2 . For the 300-day mission, 99. 9978% is dissipated as heat to the LH 2 . 
The results are similar for the other five tanks. Therefore, it is assumed 
that all of the energy input to the electric pump ends up heating the LH 2 and 
causing "boiloff. " 
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APPENDIX G 


V 


TANKAGE INSULATION OPTIMIZATION AND TVS 
HEAT EXCHANGER DESIGN 


For a given tank and mission, assuming constant tank pressure, the 
optimum insulation system can be determined by minimizing the combined 
weight of the multilayer insulation and the LH 2 "boiloff" due to external heat 
leak. First, the insulation effective conductivity must be determined. 

Based on extensive experimental work by MDAC (ref. 28), the performance 
of a typical high performance insulation, namely double- aluminized mylar 
with B4A dacron-net spacer, is characterized by the equation: 


-> r _ o / ( T w + T ) 

K , = 3. 007 x lO"^ 3 N ' — — 

eff 2 


8 . 333 x 10 ' 2 <t (T 2 + T 2 ) (T„ + T ) N 

ri O ri U 


(N-l) (4- + T~ 
1 2 


1) N 


(G-l) 


in English engineering units, where N is the layer density (assumed for our 
study at 100 layer-pairs/inch), N is the number of layer-pairs, er is the 
Stefan- Boltzmann constant, 0. 1714 x 10 - ^ Btu/hr-ft^- °R"*, Tpj = 460 °R, 


T 


q = 40°R, e 1 = e 2 = 0.021. Therefore, for N >100 


K .. = 2. 1 Joule/m-sec- °K 
eff 

(1.351 x 10 " 5 Btu/hr-ft- °R) 


(G-2) 


This value must be further degraded by 50% to account for heat leak 
through j oints, fasteners, perforations, etc. The weight of this insulation 
in kg is 

W = 0. 0122 ^ 5 - x 3, 937 S x « (A ) (m 2 ) 

n ■ \ 2 m T 

(layer-pair) m 

+ 0.229 (^) A x (m 2 ) (G-3) 

m 


where ? is the insulation thickness, (m) and A^ the insulation area (m ) 
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Thus, the total weight (boiloff plus insulation weight) is: 


W 


1. 5 

q 

+ 


( 2 . 1 ) 
e 


(255 - 22. 2) 


Aj t (sec) 


449, 000 Joules/kg 
t 


+ 48fi Aj. + 0. 229 A ] . 


HEAT SHORT 
449,000 Joules /kg 


(G-4) 


Differentiating with respect to C and equating to zero gives: 


or 


W 

dfi 


-3 "^T ^ 

-1. 64 x 10 5 — y — + 48 A t = 0 
£2 I 

(G-5) 

n 2 1. 64 x 10 -3 , 

“ 48 t ' 

(G-6) 

8 _ 5. 84 x 10“ 3 1/2 

(G-7) 


i 

Solving equation (G-7) for the insulation thickne'ss (m) in terms of 
mission time (sec), t, the boiloff flowrate is: * 


W 


1.64x10 Aj ^HEAT SHORT' 
e + 449,000 


(G-8) 


This boiloff, together with that caused by pump/motor input power, 
gives the total vented flowrate through the TVS heat exchanger. By examining 
the heat transfer processes which occur, it will be possible to define the 
required heat transfer area and heat exchanger size. 

Assuming the configuration shown in Figure 104 with a heat- exchanger 
coil bonded or brazed to the outside of the standpipe (so as not to affect 
flow or pressure-drop in the standpipe), the required length of coil and heat 
exchanger area will be defined by the overall heat-transfer coefficient. 


The flow through the standpipe is being pumped by the TVS pump at 1% 
or 0. 1% tank volume /minute; thus, the heat transfer coefficient on the inside 
of the standpipe, h^, is governed by forced convection. Forced convection 
in a circular duct is described by the Dittus- Boelter equation: 



(G-9) 
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STANDPIPE 


Figure 104. TVS Heat Exchanger Nomenclature 

• o 

where Q is the TVS flowrate in mV /sec, D g is the standpipe diameter in m, 
and h^ is in Joule/m^- sec- °K with LH 2 properties evaluated at 29. 2° K. 
Inside the heat exchanger tubing the heat transfer coefficient, h^, is also 
governed by forced convection because of the vent flow, and the same 
equation applies: 


h 2 P TUBE 

K 


= 0.023 


4W 


0 . 8 




Pr 


0. 33 


TUBE 


(G- 10) 


where W is the vent flow in kg/sec, DtUBE i- s the heat exchanger tube 
diameter in m and is in Joules /m.2 - sec- °K with LH 2 properties evaluated 
at a film temperature of 21. 1° K. 

Referring to Figure 104, assume that the tubing is spaced 0.318 cm 
(1/8 inch) apart, and only the half of the heat exchanger tubing in contact 
with the TVS flow is used for heat transfer to the flow in the tubing. The 
outside of the heat exchanger will be insulated with foam insulation to pre- 
vent condensation from the tank ullage which could use up all of the cooling 
capacity of the vent fluid. The foam insulation thickness will be sized to 
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limit the condensing heat transfer from the ullage to 10% of that from the 
TVS flow. The heat flux required is that which will boil the vent flow, and 
equating this heat flux to the overall heat transfer coefficient, including the 
standpipe wall and tube resistance, K/X, gives: 


(1. 1) W (449. 000 Joules/kg) = 


L (T H z " t hex ) 


_J , X 1 

h : D 2 KD 2 h 2 jr(D 2 -2t)/2 


(G- 1 1) 


where 


X 


(t 4 D 2 /2) D 2 

d: 


tt(D 2 - 2t') 

8 


or the heat exchanger length, L, is. 


L 


W (449. 000) (1. 1) 
( t h 2 " t hex) 


h l°2 


(t + D 2 /2 - 7r(D 2 - 2t») /8D 

kTx 


+ 


2 

h 2 7r (D 2 - 2t) 


(G- 12) 


To minimize L, the smallest practical value of (tube diameter) 
should be used; however, the tube must be large enough to give a maximum 
pressure drop of about 1. 38 N/cm^ (2 psi) at the vent flow W. The value of 
1. 38 N/cm^ (2 psi) was chosen because it is assumed that the vent flow is 
expanded to Tjqjr;x = 17. 2 K (31° R), (Th 2 = 25. 2° K (45. 4° R)) at 3. 45 N/cm^ 
(5 psia). With a pressure drop in the tube of onlv 1. 38 N/cm.2 (2 psi), the 
vent flow will exit the system at about2- 07 N/crm (3 psia), which is enough 
above the LH 2 triple-point pressure of 0. 69 N/cm^ (1 psia) to preclude 
premature freezing of the vent flow upon expansion to space. 


The pressure drop in the tube is 


AP 



W 


2 pg c (f d 2 ¥ 


(G-13) 
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The friction factor, f, can be evaluated for turbulent flow in a smooth 
tube from the Blasius correlation 


0. 316 
r0.25 


(G- 14) 


where the Reynolds number, R, is that evaluated for the vent flow in 
equation (G-10) 

The foam insulation thickness, Xp, can be determined from the heat 
flux criteria described above, and the equation: 


0. 1 (W) 449, 000 = 


L (t t “ t hex> 


X, 


(G-15) 


K f (D 2 + D 3 ) h 2 tt(D 2 - 2t') 


The heat exchanger weight is found from: 


W 


HEX 


D_ 


n D_ 


L D 2 n t' + ' 


HEX 


(G- 16) 


The heat exchanger length, tube diameter, weight, foam thickness and 
foam weight for each of the six tanks are shown in Table XIV. Since the 
thickness for the foam is too small for practical fabrication, 1.27-cm (1 72- 
inch) of foam was arbitrarily used. The foam weight shown is for the 
1.27 cm (1 /2-inch) thickness. 


197 



TABLE XIV. - TVS HEAT EXCHANGER DESIGN PARAMETERS 


Tank ID 

Tube 

Tube wall 


Heat exchanger 

Foam 

Foam 


diameter. 

thickness. 

Tube length. 

weight. 

thickness, a 

weight k _ 


cm (in. ) 

cm (in. ) 

m (ft) 

kg (lb) 

cm (in ) 

kg (lb) 

5, 000/4 

0. 635 (0. 25) 

0.0508 (0.02) 

6. 9 (22. 6) 

0. 577 (1. 275) 

0. 104 (0. 041) 

0.0385 (0.085) 

500/4 

0. 318 (0. 125) 

0.0305 (0.012) 

2. 59 (8. 5) 

0. 066 (0. 145) 

0. 094 (0. 037) 

0. 0109 (0. 024) 

500/2 

0. 318 (0. 125) 

0. 0305 (0. 012) 

2.44 (8.0) 

0. 063 (0. 140) 

0. 107 (0. 042) 

0. 0104 (0. 023) 

500/1 

0. 318 (0. 125) 

0.0305 (0. 012) 

2. 38 (7. 8) 

0. 061 (0. 135) 

0. 114 (0. 045) 

0.01 (0.022) 

50/2 

0. 16 (0. 063) 

0.0152 (0. 006) 

0. 915 (3. 0) 

0. 00425 (0. 0094) 

0. 140 (0. 055) 

0. 00385 (0.0085) 

50/1 

0. 16 (0. 063) 

0.0152 (0. 006) 

0. 885 (2. 9) 

0.0040 (0.0089) 

0. 145 (0. 057) 

0.0037 (0. 0082) 

a Based on 10% of design 
^Based on 1 . 27 - cm (0.5 

heat flux through foam 
in.) thickness 
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